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INTRODUCTION 


I.  SOLID  STATE  DEVICE  RESEARCH 

Photovoltaic  response  has  been  observed  in  p-n  junctions  of  Fbj_xSn  Te  at  wavelengths  up  to 
1 1  \i  at  77°K  and  14  \i  at  12  °K.  Noise  measurements  at  12  °K  give  a  noise  equivalent  power  for 

-12  1/2  Q  1/2 

a  Fb^  g^Sn^  j^Te  diode  of  about  5  x  10  W  sec  '  and  a  D*  of  8  *  10  cm/W  sec  '  at  ll|i. 
The  responsivity  cutoff  of  a  Pb^  g^Sn^  ^Te  diode,  taken  at  50  percent  of  the  peak  respon- 
sivity,  is  14 ^  at  12 °K  and  1 1  jjl  at  77°K.  These  results  indicate  that  alloys  in  this  series  have 
considerable  potential  for  infrared  detection  throughout  the  8-  to  14 -|jl  atmospheric  window  and 
well  beyond. 

Infrared  laser  emission  has  been  observed  at  12 °K  and  77°K  from  diodes  of  Pb.  Sn  Se  and 

1-y  y 

Pb  Sn  Te.  A  Pb  Sn  Se  diode  with  a  nominal  Sn  composition  of. 5  percent  produced  laser 

emission  at  1 2 .4  jjl  at  12°K  and  at  10.2  fa  at  77°K.  Spontaneous  emission  from  different 

Pb,  Sn  Te  diodes  at  12° K  was  observed  at  9.4,  12.7,  13.5  and  13.7 jjl  The  third  diode  in  this 
1  -x  x 

series  also  produced  laser  emission  at  10.6 [i  at  77° K. 

The  Interface-alloy  technique  has  been  used  to  produce  heterojunctions  between  GaAs  and 
InSb.  Despite  the  relatively  large  14  percent  lattice  mismatch  between  the  semiconductors, 
these  heterojunctions  are  single-crystal,  with  an  interdiffusion  of  the  four  elements  of  less 
than  2  \±.  Photocurrent,  current-voltage  and  capacitance-voltage  measurements  are  explained 
by  a  model  for  the  heterojunction  band  structure  employing  a  linearly-graded  energy  gap  of 

o 

the  order  of  60A  long  joining  the  GaAs  to  the  InSb. 

A  400  keV  Van  de  Graaff  positive  ion  accelerator  has  been  used  to  implant  phosphorus  ions  in 

germanium.  Experimental  evidence  indicates  that  the  implanted  ions  have  produced  n-type 

regions  in  bulk  p-type  material.  Implantation  for  one  hour  with  an  incident  phosphorus  beam 

current  of  0.4  pA  resulted  in  a  phosphorus  concentration  at  the  sample  surface  of  the  order  of 
20  -3 

1  x  10  cm  ,  and  an  n-type  region  approximately  2  pdeep. 

By  taking  the  negative  mobility  as  an  adjustable  parameter,  we  have  obtained  good  quantitative 
agreement  between  the  measured  small -signal  admittance  parameters  of  GaAs  negative  con¬ 
ductance  amplifiers  and  those  calculated  from  the  three-slope  piecewise-linear  model  for 
velocity  vs  electric  field.  The  same  set  of  parameter  values  also  provides  both  a  reasonable 
fit  to  the  measured  current-voltage  characteristics  and  a  correct  prediction  of  the  frequency 
and  threshold  bias  for  oscillation. 

Very  pure  epitaxial  layers  of  GaAs  have  been  grown  on  Cr-doped  semi-insulating  GaAs  sub¬ 
strates .  The  best  layers  have  electron  concentrations  in  the  high  lO1^  to  low  lO^cm  ^ 

range,  both  at  room  temperature  and  at  77° K.  The  electron  mobilities  for  these  layers  are 
2 

typically  8000cm  /Vsec  at  room  temperature,  and  increase  at  77° K  to  values  as  high  as 
2 

77,000cm  /Vsec.  Gunn  effect  oscillators  made  from  similar  films  deposited  on  heavily 
doped  substrates  have  yielded  power  conversion  efficiencies  as  high  as  9  percent  from  de¬ 
vices  delivering  a  peak  output  power  of  1.75  watts  at  9.3 GHz.  This  represents  a  substantial 
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improvement  over  previously  reported  results  of  other  workers  and  demonstrates  that  Gunn 
oscillators  with  relatively  high  efficiencies  can  be  fabricated  using  epitaxially  grown  GaAs. 

11.  OPTICAL  TECHNIQUES  AND  DEVICES 

Heterodyne  sensitivity  at  10.6  [xm  has  been  measured  using  a  Cu-doped  germanium  photocon¬ 
ductor  at  liquid  helium  temperature.  The  measured  minimum  detectable  power  is  within  10 dB 

-20 

of  the  theoretical  limit,  h^B  =  2  x  10  W/Hz  bandwidth.  The  signal  frequency  was  60kHz 
and  the  noise  bandwidth  approximately  300  kHz. 

A  sealed-off  CO^  laser  has  been  operated  for  approximately  100  hours  without  serious  loss  in 
output  power.  Although  the  gas  fill  was  CO^,  ^  and  He,  spectra  of  the  visible  spontaneous 
emission  from  the  side  of  the  tube  indicate  the  presence  of  CO,  CN,  and  N^-  The  mechanism 
of  formation  of  the  additional  constituents  is  being  investigated  to  determine  their  relation  to 
the  operating  life  limitations  of  sealed-off  systems. 


111.  MATERIALS  RESEARCH 

The  crystal  structure  of  Nd^O^Te  has  been  determined  by  computer  analysis  of  x-ray  powder 
diffraction  data.  The  atomic  arrangement  is  derived  from  the  K^NiF^  structure,  in  which  the 
rare  earth  atoms  exhibit  an  unusual  nine-fold  coordination. 

In  order  to  determine  the  cause  of  the  anomaly  observed  in  the  x-ray  scattering  of  ZnSe, 
neutron  diffraction  measurements  have  been  made  on  this  compound  at  room  temperature  and 
liquid  helium  temperature.  The  results  obtained  so  far  are  inconclusive  because  of  difficul¬ 
ties  in  evaluating  background  intensity  corrections. 

The  phase  diagram  of  Ag^Te  at  high  pressure  has  been  investigated  by  means  of  electrical 
resistivity  and  x-ray  diffraction  measurements.  Evidence  has  been  obtained  for  three  high- 
pressure  phases,  two  of  which  exist  at  room  temperature. 

Hysteresis  of  the  first-order  B8j  ^  B31  magnetic  transition  in  MnAs  has  been  studied  by 
measuring  the  resistivity  and  magnetic  susceptibility  as  functions  of  pressure  and  tempera¬ 
ture.  The  experimental  results,  together  with  earlier  data  for  MnASj^P^,  show  that  the 
transition  is  due  to  a  large  exchange  striction  in  the  basal  planes  together  with  a  volume- 
dependent  Weiss  molecular  field  and  manganese  atomic  moment. 

The  retrograde  solubility  of  Pb  in  PbS  has  been  investigated  by  means  of  Hall  coefficient 
measurements  on  n-type  samples  prepared  by  annealing  and  quenching  single  crystals  con¬ 
taining  excess  Pb.  The  carrier  concentration  decreases  with  decreasing  annealing  tempera  - 

1  Q  _  Q  lO  _0 

ture  from  2.4  x  10  cm  for  912°C  to  1.3  x  10  cm  for  350°C. 

Measurements  of  the  Hall  coefficient  vs  temperature  for  S-doped  GaSb  have  shown  that  the 
lowest  donor  level  of  sulfur  lies  about  0. 075 eV  below  the  lowest  conduction  band  at  atmos¬ 
pheric  pressure.  The  pressure  dependence  of  resistivity  at  room  temperature  shows  that 
this  donor  level  is  associated  primarily  with  the  conduction  band  minima  at  the  (100)  zone 
faces . 
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The  magnetoresistance  and  Hall  coefficient  of  high-purity  and  V-doped  T^O^  single  crystals 
have  been  measured  at  magnetic  fields  up  to  218  kG.  The  high  magnetoresistance  of  the  high- 
purity  samples  at  4.2°K  is  inconsistent  with  a  conventional  one-band  carrier  transport  mech¬ 
anism.  The  p-type  conductivity  and  negative  magnetoresistance  of  V-doped  suggest 

that  part  of  the  V  is  present  in  the  quadrivalent  state. 

Efficient  pulsed  laser  operation  has  been  obtained  in  Nd-doped  YVO^  single  crystals .  The 
threshold  is  only  about  1  joule.  This  low  threshold  is  due  to  the  reduced  Stark  splitting  of 
Nd4^  ions  in  YVO^  and  to  substantial  transfer  of  absorbed  energy  from  the  lattice  to  the  Nd+'^ 
ions . 

A  method  employing  PbTe  and  SnTe  as  standards  has  been  developed  for  determining  the 
presence  of  Pb  and  Sn  inPbj_xSnxTe  alloys  by  electron  microprobe  analysis.  The  results  are 
in  good  agreement  with  those  obtained  by  wet  chemical  techniques. 


IV.  PHYSICS  OF  SOLIDS 

The  optical  study  of  metals  and  semiconductors  by  means  of  sensitive  differential  techniques 
is  continuing.  One  of  these  methods,  the  current  pulse  modulation  technique,  has  been  ap¬ 
plied  to  the  study  of  gold  and  nickel  in  the  energy  region  from  2  to  10 eV;  much  well  defined 
structure  not  observed  in  the  standard  reflectivity  has  been  found.  This  same  modulation 
technique  is  also  being  used  in  InSb  to  study  the  magneto-oscillations  associated  with  transi¬ 
tions  across  the  direct  energy  band  gap,  and  also  from  the  spin-orbit  split -off  valence  band  to 
the  k  =  0  conduction  band. 

Oscillations  in  the  infrared  magnetoreflectivity  of  arsenic  have  been  observed  for  the  first 
time.  Two  distinct  sets  of  interband  transitions,  taking  place  at  different  points  in  the 
Brillouin  zone,  have  been  seen. 

The  optical  study  of  the  effects  of  calibrated  uniaxial  stress  on  infrared  transitions  inS-doped 
silicon  continues.  A  value  for  the  pure  shear  deformation  potential  of  the  conduction  band, 
E2  =  7.9 eV,  is  obtained  from  the  behavior  of  the  p-like  transitions.  The  application  of  stress 
helps  distinguish  between  possible  symmetries  of  observed  centers,  andaids  in  identification. 

While  studying  polaron  effects  in  InSb,  a  discontinuity  in  the  functional  dependence  of  the 
polaron  cyclotron  energy  on  magnetic  field  was  observed  near  the  LO  phonon  energy;  this 
phenomenon  had  been  predicted  theoretically.  In  addition,  absorption  lines  associated  with 
donor  states  were  found  to  exhibit  polaron  effects. 

Experiments,  in  collaboration  with  Professor  W.  Paul's  group  at  Harvard,  to  establish  the 
band  ordering  in  the  mercury  chalcogenides  have  now  been  extended  to  studies  of  the  pres¬ 
sure  dependence  of  the  absorption  edge.  These  new  experiments  suggest  the  same  relative 
energy  level  orderings  at  k  =  0  for  the  and  bands  in  HgSe  and  HgTe,  which  is  inconsist¬ 
ent  with  the  previous  conclusion  derived  from  the  pressure  dependence  of  the  thermoelectric 
power.  This  contradiction  appears  to  have  been  resolved  by  Hall  measurements  under  pres¬ 
sure  which,  although  incomplete,  indicate  that  at  least  some  of  the  HgSe  samples  were  not 
extrinsic,  thereby  invalidating  the  thermoelectric  power  interpretation. 
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An  analysis  is  being  made  of  the  beating  effects  observed  by  Whitsett  in  Shubnikov-de  Haas 
measurements  in  n-type  HgSe.  Using  the  theoretical  energy  momentum  expressions  derived 
by  Kane  for  zinc  blende  materials  with  small  energy  gaps,  an  attempt  is  being  made  to  con¬ 
struct  the  Fermi  surface  for  this  semiconductor. 

The  conductivity  tensor  in  bismuth  has  been  calculated,  assuming  an  anisotropic  effective 
mass  for  holes  and  the  two-band  model  for  the  electrons.  A  detailed  lineshape  calculation  for 
the  magnetoreflection  is  now  under  way. 

The  electronic  energy  bands  in  silicon  have  been  calculated  by  using  the  Fourier  expansion 
technique  in  which  the  band  parameters  are  determined  from  the  experimental  values  for  the 
energy  gaps  and  effective  masses.  This  band  structure  has  allowed  the  calculation  of  a  fre¬ 
quency  dependent  dielectric  constant  which  is  in  reasonably  good  agreement  with  experiment. 

Calculations  using  high  temperature  expansion  techniques  for  an  ensemble  of  spins  interacting 
with  a  Heisenberg  exchange  interaction  are  continuing.  Thermodynamic  functions ,  such  as 
the  spin  correlation  function  and  zero  field  susceptibility,  are  being  treated  for  any  range  of 
exchange  interaction  and  any  lattice.  Higher  order  coefficients  of  the  high  temperature 
series  for  the  zero  field  susceptibility  have  been  incorporated  in  the  arguments  for  phase 
transitions  in  two-dimensional  Heisenberg  models. 

The  study,  in  the  time  dependent  Hartree-Fock  approximation,  of  spin  sound  (or  spin  density 
waves)  in  unmagnetized  or  weakly  magnetized  Fermi  gases  has  been  concluded;  the  results 
will  appear  in  the  Physical  Review.  As  a  sequel  to  this  work,  the  properties  of  ordinary  den¬ 
sity  waves  have  been  studied  in  the  same  approximation.  In  addition  to  the  usual  plasmon 
iispersion  relation,  another  branch  of  the  density  oscillation  spectrum  is  discussed. 

In  studying  the  influence  of  the  electron-electron  interaction  on  the  transport  properties  of  the 
conduction  electrons  in  a  metal,  a  model  of  an  interacting  electron  gas  driven  by  an  external 
electromagnetic  field  and  scattered  by  randomly  distributed  impurities  has  been  considered; 
the  electron-electron  interaction  is  treated  in  the  random  phase  approximation.  So  far  this 
theory  has  been  applied  to  the  special  case  of  a  homogeneous  electron  gas  in  the  Born 
approximation . 

For  the  case  of  a  Landau  quasi-particle  distribution  function  in  the  presence  of  random  im¬ 
purities,  the  transport  theory  is  incomplete  since  it  is  only  valid  for  go  ^  0  or  c or  «  1;  for 
c or  »  1,  the  case  of  experimental  interest,  the  result  is  incorrect.  The  possibility  of  a  new 
approach,  valid  for  an  arbitrary  cj  is  now  being  investigated. 

Further  results  on  stimulated  scattering  experiments  in  quartz  at  2.1°K  have  lent  additional 
support  to  the  contention  that  the  anomalous  Brillouin  shifts  reported  by  Krivokhizha  et  al . 
were  misinterpreted  and  are  actually  due  to  intense  Raman  radiation  which  sets  in  at  low 
temperatures.  This  phase  of  the  experiment  is  now  completed;  results  have  been  submitted 
to  Science. 

A  study  of  the  effect  of  orientational  correlation  and  molecular  fields  in  second  harmonic  light 
scattering  from  liquids  and  gases,  with  application  to  temperature  dependence  and  depolariza¬ 
tion,  has  been  concluded.  The  results  are  being  prepared  for  publication. 
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Work  on  infrared  Raman  scattering  in  solids  continues.  Raman  spectra  have  been  observed 
in  GaS,  GaSe  and  trigonal  selenium. 

The  analysis  of  a  maser  model  of  interstellar  OH  microwave  emission  has  been  published  in 
Physical  Review  Letters.  Further  calculations  of  self-consistent  populations  under  direc¬ 
tional  UV  pumping  and  the  effect  of  microwave  saturation  on  mode  suppression  are  now  being 
carried  out. 
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I.  SOLID  STATE  DEVICE  RESEARCH 


A.  PHOTOVOLTAIC  EFFECT  IN  Pb^Si^Te  DIODES 

We  have  observed  photovoltaic  response  in  p-n  junctions  of  Pb^^Sn^Te  at  wavelengths  up 
to  11  p  at  77  °K  and  14  p  at  12°K.  These  results,  together  with  previous  photoluminescence  data 
and  the  proposed  model  for  the  band  structure  of  Pb^  ^Sn^Te  (Ref.  1),  indicate  that  these  alloys 
have  considerable  potential  for  efficient  infrared  detection  throughout  the  8  to  14  p  atmospheric 
window  and  well  beyond. 

The  results  reported  here  were  obtained  with  p-n  junctions  which  formed  in  the  vapor  dep¬ 
osition  process  used  in  growing  the  alloys.  Subsequently,  junctions  were  also  prepared  by  a 

2 

diffusion  technique  previously  used  in  fabricating  PbTe  diode  lasers.  The  alloy  crystals  had 
carrier  concentrations  in  the  range  of  2  to  5  x  10^7  cm  ^  and  mobilities  in  excess  of  10  cm"/ 

V  sec  at  77  °K.  Figures  I-l(a)  and  -1(b)  show  current-voltage  characteristics  of  a  Pb^  ^Sn^  ^Te 
and  of  a  Pb^  ^Sn^  ^Te  diode  respectively,  at  12°K.  These  diodes  were  in  the  form  of  rec¬ 
tangular  parallelepipeds  with  cleaved  surfaces  normal  to  the  plane  of  the  junction.  The  samples 
were  mounted  on  a  copper  heat  sink  with  the  p-type  layer  attached  to  the  copper.  The  thin 
n-type  layer  (<20p)  was  exposed  to  the  signal  light.  Liquid  solder  was  used  for  making  contacts 
in  order  to  avoid  heating  the  sample.  In  these  early  diodes  no  effort  was  made  to  treat  the 
surface  by  etching.  The  reverse  current  in  Figs.  I-l(a)  and  -1(b)  is  thought  to  be  mostly  due 
to  surface  leakage.  At  77 °K  the  reverse  current  in  these  diodes  was  somewhat  higher  than 
at  12°K. 


(a)  PbQ,  91  Sng.  09^e  diode  with  area  (b)  Pbg  86^n0.  14"^e  diade  with  area 

af  0.  33  X  0.  20  mm2.  0f  q  28  X  0.  58  mm2. 


Fig.  1-1.  Current  val  tage  characteristics  at  12°K. 
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Section  I 


Fig.  1-2.  Responsivity  spectro  of  Pb(j  ^^SnQ^  ^Te 
diode  of  Fig.  1-1  (b)  at  12°K  and  ot  77° K. 


Fig.  1-3.  Responsivity  spectra  of  Pb^  o'5^nn  17T© 
diode  ot  12°K  and  ot  77°K.  *  * 


x 


Fig.  1-4.  Response  of  Pbo. 83Sn0.  17Te 
diode  ot  77°K  to  o  Q-switched  CO^ 
loser  pulse. 


Fig.  1-5.  Energy  gop  of  Pb-|_xSnxTe  os  o  function  of  mole 
froction  of  SnTe,  x.  The  doto  points  ore  obtoined  from 
the  cutoff  of  the  photovoltoic  responsivity  (this  work)ond 
from  luminescence  spectro. 
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Section  I 


Figure  1-2  shows  the  spectral  dependence  of  the  open-circuit  responsivity  obtained  with  the 
PbQ  86^n0  14^e  d^oc^e  zero  bias.  The  responsivity  at  12°K  has  a  value  of  57  V/W  at  11  p  and 
at  77°K  a  value  of  1.7  V/W  at  9  p.  The  enhancement  observed  in  Fig.  1-2  at  both  temperatures 
just  prior  to  the  long-wavelength  cutoff  can  be  attributed  to  an  increase  of  carrier  density  at 
the  junction  as  the  radiation  begins  to  penetrate  the  n-type  layer  in  the  vicinity  of  the  bandgap 
absorption  edge  of  the  crystal.  For  wavelengths  shorter  than  5  p,  the  responsivity  at  both  tem¬ 
peratures  has  a  slope  slightly  greater  than  unity,  the  slope  expected  for  an  ideal  photon  detector. 

The  number  of  carriers  reaching  the  p-n  junction  per  photon  incident  on  the  surface  of  the 
n-type  layer  has  been  estimated  to  be  0.1  at  12°K  for  the  lip  wavelength.  This  estimate  was 
based  on  a  measurement  of  the  short-circuit  photocurrent.  The  fractional  loss  of  0.9  consists 
primarily  of  a  loss  of  about  0.5  due  to  a  reflection  at  the  diode  surface,  and  of  losses  due  to 
recombination  at  the  surface  and  in  the  bulk  of  the  n-type  layer.  No  attempts  have  been  made 
to  reduce  these  factors  in  the  present  devices. 

The  speed  of  the  photovoltaic  response  was  measured  by  means  of  light  pulses  from  a  GaAs 
diode  laser,  which  had  a  risetime  of  less  than  1  nsec.  The  time  constant  of  the  photovoltage 
pulse  observed  by  connecting  a  properly  terminated  50-ohm  line  directly  to  the  detector  diode 
was  about  200 nsec  at  12°K  and  about  20 nsec  at  77 °K.  The  200  nsec  value  provides  an  estimate 
of  the  carrier  lifetime  at  12°K.  The  20 nsec  value  is  an  upper  limit  to  the  77 °K  lifetime,  since 
at  20  nsec  the  capacitance  of  the  diode  could  be  a  limiting  factor  to  the  response.  The  lower 
lifetime  at  77 °K  may  partly  account  for  the  reduced  responsivity  at  this  temperature,  as  shown 
in  Fig.  1-2,  because  the  fraction  of  generated  carriers  which  *  ch  the  junction  is  directly  re¬ 
lated  to  the  diffusion  length. 

The  rms  noise  voltage  at  12 °K  was  found  to  be  approximately  proportional  to  the  reciprocal 

of  the  frequency  between  500  and  3000  Hz.  At  the  900  Hz  chopping  frequency,  where  the  above 

-10  1/2 

responsivity  measurements  were  made,  a  noise  voltage  of  3  X  10  V  sec  '  was  determined 

-12  1/2 

at  12°K.  This  corresponds  to  a  noise  equivalent  power  of  about  5  X  10  W  sec  '  and  a  D* 

9  ,  1/2 

of  8  x  10  cm/W  sec  '  at  the  lip  wavelength. 

Figure  1-3  shows  the  spectral  dependence  of  the  open-circuit  responsivity  of  a  PbQ  g^SnQ  ^Te 

diode  at  zero  bias.  The  responsivity  and  D‘  for  this  diode  are  about  the  same  or  slightly  larger 

than  for  the  previous  diode.  The  responsivity  cutoff,  taken  at  50  percent  of  the  peak  responsivity 

is  14  p  at  12°K  and  11  p  at  77° K.  This  means  that  at  12°K  this  diode  responds  throughout  the 

atmospheric  window  and  at  77 °K  it  will  detect  the  CO^  laser  radiation  at  10.6  p.  The  time 

response  of  the  PbQ  ^Sn^  ^Te  detector  at  77 °K  to  a  Q-switched  C02  laser  pulse  is  shown  in 

Fig.  1-4.  Figure  1-5  shows  the  variation  of  the  energy  gap  with  alloy  composition  at  both  12°K 

12  3 

and  77 °K  as  determined  by  luminescence  data  *  *  and  by  the  photovoltage  measurements.  The 

temperature  coefficient  of  the  energy  gap  of  the  alloys  between  12°K  and  77 °K  agrees  closely 

with  the  temperature  coefficient  in  PbTe,  and  the  energy  gap  at  77 °K  decreases  with  increased 

Sn  concentration  with  a  slope  nearly  equal  to  the  slope  at  12 °K.  These  observations  lend  further 

1 

support  to  the  model  for  the  band  structure  of  the  alloys. 

I.  Melngailis  T.C.  Harman 

A.  R.  Calawa  J.O.  Dimmock 

J.  F.  Butler 
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B.  DIODE  LASERS  OF  Pb1  Sri  Se  AND  Pb1  Sn  Te 

i  y  y  i “X  x 

We  should  like  to  report  the  observation  of  infrared  laser  emission  at  12°K  and  77 °K  from 
diodes  of  Pb  Sn  Se  and  Pb  Sn  Te.  Optically  pumped  infrared  laser  emission  in  Pb  Sn  Te 

*-y  y  lxx|  45  1  "x  x 

at  12°K  has  recently  been  observed.  From  those  and  other  experiments,  '  it  was  concluded1 
that  the  energy  gap  of  Pb^  ^Sn^Te  at  12°K  decreases  smoothly  with  increasing  Sn  concentration, 
becoming  zero  at  x  =  0.35.  Our  current  results  show  that  the  energy  gap  of  Pb^  Sn^Se  also 
decreases  with  a  decreasing  Pb:Sn  ratio,t  and  indicate  that  the  value  of  y  for  zero  energy  gap 
at  12°K  lies  within  the  range  0.11  ^  y  ^  0.14.  This  work  demonstrates  that  diode  laser  action 
over  a  wide  range  in  the  infrared  can  be  attained  with  the  (Pb,  Sn)  salts. 

Pb^  ySn^Se  crystals  were  vapor-grown  at  825°C,  using  a  10-gm  powdered  source  ingot  of 
composition  y  =  0.05,  and  air-cooled  from  the  growth  temperature.  The  source  ingot  was  pre¬ 
pared  from  a  100-gm  stoichiometric  mixture  of  the  components  that  had  been  heated  above  the 
melting  point  and  water-quenched.  Vapor-grown  crystals  and  the  quenched  ingot  were  n-type. 

A  diffused  p-type  skin  was  produced  by  heating  the  crystals  for  one  hour  at  400°C  with  a  vapor 
source  of  Se-saturated  PbSe.  Crystal  growth,  ingot  preparation  and  diffusion  runs  were  carried 
out  in  evacuated,  sealed  quartz  ampules.  Electron  beam  microprobe  analysis  yielded  a  value 
of  y  =  0.032  on  the  diffused  surface  and  y  -  0.044  in  the  crystal  bulk.  In  the  p-n  junction  region, 
y  is  expected  to  have  an  intermediate  value.  Diodes  in  the  form  of  rectangular  parallelepipeds 
were  made  by  cleaving  along  { 1 00 ]  surfaces,  and  were  mounted  in  low  inductance  packages, 
using  a  low  melting  temperature  Hg-In-Tl  solder. 

Pbj  xSnxTe  crystals  for  several  values  of  x  were  vapor  grown  by  a  similar  method.  Growth 
temperatures  were  between  700°C  and  825°C.  For  nominal  values  of  x  between  0.16  and  0.20 
p-type  vapor  grown  crystals  were  obtained  using  metal  saturated  source  powders.  An  n-type 
skin  formed  during  the  cooling  process.  For  x  0.10,  n-type  Pb^_xSnxTe  was  vapor  grown 
from  metal  saturated  source  powders  and  a  p-skin  was  formed  subsequently  by  diffusion  using 
a  Te-saturated  PbTe  source.  Diode  lasers  have  also  been  fabricated  from  Bridgman-grown 
Pb.  ,Sn  Te.  The  Pb  Sn  Te  diodes  were  shaped  and  packaged  in  the  same  manner  as  were 
the  Pb,  Sn  Sc  devices. 

*-y  y 

For  spectral  measurements,  diodes  were  mounted  on  the  cold  finger  of  a  liquid  helium 

dewar  equipped  with  KRS-5  windows.  The  diodes  were  pulse  excited,  although  CW  operation 

at  12°K  was  possible  in  both  Pb,  Sn  Se  and  Pb,  Sn  Te  devices.  Spectra  were  measured 

1  -y  y  1  -x  x 

using  a  NaCl  prism  monochromator,  a  copper-activated  germanium  detector  and  a  lock-in 
amplifier  system. 

Figure  1-6  shows  the  12°K  spectra  for  a  Pb^  ^Sn  Se  diode  at  current  levels  below  and 
above  the  threshold  for  stimulated  emission.  Threshold  current  density  at  this  temperature 
was  670  A  cm  .  The  Fabry-Pcrot  cavity  length  was  315  p.  Using  this  value  the  observed 
mode  spacing  gives  a  value  of  6.4  for  the  quantity  (n^  —  XQ  dn/dX).  At  a  temperature  of  77 °K, 
the  threshold  current  density  was  2700  A  cm  ^  and  laser  action  occurred  at  10.2  p.  The  emis¬ 
sion  data  of  Fig.  1-6  indicate  an  energy  gap  at  12°K  of  O.lOOeV  in  the  p-n  junction  region  of  the 
Pb^  ^Sn^Se  device,  while  the  energy  gap  of  pure  PbSe  at  this  temperature  is  0.145eV.^  This 


t  A.  J.  Strauss  has  previously  obtained  this  result  from  optical  absorption  measurements  on  thin  films  of  Pbj  ^SnySe  for  0  <  y  <  0.20. 
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Fig.  1-6. 


Emission  spectra  far  Pb^  ^Sn^Se  diade  laser  at  approximately  12°K. 


dependence  of  energy  gap  on  Sn  concentration  for  Pb^  ySn^Se  is  similar  to  that  for  Pb^  ^Sn^Te. 
Linear  extrapolations  using  the  energy  gap  of  0.145eV  at  y  =  0  and  0.100  eV  at  the  various  y- 
values  found  by  electron  beam  microprobe  measurements  on  the  laser  diode  indicate  a  zero 
energy  gap  at  12  °K  in  Pb^  ^Sn^Se  for  some  value  of  y  between  0.11  and  0.14. 

Spontaneous  and  coherent  emission  spectra  for  a  Pbo#83SnQ  ^Te  diode  at  12°K  are  given 
in  Fig.  1-7.  The  threshold  current  density  is  600  A  cm"  .  The  cavity  end  faces  were  220  p 
apart.  Using  this  value  the  mode  spacing  gives  a  value  of  6.9  for  (nQ  —  X^dn/dX).  At  77 °K 
laser  action  at  10.6  p  was  observed  with  a  threshold  current  density  of  3,  000  A  cm  ^ .  Diode 
laser  action  at  12°K  for  other  values  of  x  in  the  Pb,  Sn  Tc  system  was  observed  at  13.7, 

12.7  and  9.4  p.  The  dependence  of  emission  wavelength  on  Sn  concentration  in  the  Pb^  ^Sn^Te 
diodes  is  in  substantial  agreement  with  Fig.  3  of  Ref.  1. 


J.  F.  Butler 
A.  R.  Calawa 
T.  C .  Harman 


C.  THE  GaAs-InSb  GRADED-GAP  HETEROJUNCTION 

A  paper  entitled  "The  GaAs-InSb  Gradcd-Gap  Hetero junction"  has  been  submitted  to  Solid- 
State  Electronics.  The  abstract  of  the  paper  follows: 

The  interface-alloy  technique  has  been  used  to  produce  heterojunctions  between 
GaAs  and  InSb.  X-ray  and  Kossel  line  patterns  show  that,  despite  the  relatively 
large  14%  lattice  mismatch  between  the  semiconductors,  these  hetcro junctions  are 
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Fig.  1-7.  Emission  spectra  for  Pb^  ^Te  diode  laser  at  approximately  12°K. 

single-crystal.  Electron  beam  microprobe  analysis  indicates  that  the  inter¬ 
diffusion  of  the  four  elements  is  smaller  than  the  1.9  p  diameter  of  the  beam. 
Photocurrent  and  current -voltage  measurements  are  explained  by  a  model  for 
the  hetero junction  band  structure  in  which  the  salient  feature  is  a  region  of  the 

o 

order  of  60A  long  which  has  a  linearly-graded  energy  gap  joining  the  GaAs  to 
the  InSb.  In  addition,  interface  states  "fix"  the  location  of  the  GaAs  bands  at 
the  hetero  junction  interface  such  that,  at  room  temperature,  the  conduction 
band  extrapolates  to  a  value  approximately  0.93  eV  above  the  Fermi  level. 

The  heterojunction  barrier  maximum  is  smaller,  however,  because  a  portion 
of  the  GaAs  depletion  layer  potential  is  developed  across  the  graded-gap 
region.  The  photocurrent  occurs  via  hot  carriers  generated  in  the  graded-gap 

o 

region  which  traverse  this  region  (with  a  mean  free  path  of  approximately  20  A) 
to  the  heterojunction  barrier  maximum.  As  predicted  by  the  model,  incident 
monochromatic  radiation  of  energy  smaller  than  the  GaAs  bandgap  produces  a 
photocurrent  which  varies  exponentially  with  photon  energy  as  Iq  exp  [C(hy  —  E  )], 
where  C  is  a  positive  parameter  which  decreases  for  increasing  reverse  bias 
on  the  hetero  junction,  and  1^  is  the  extrapolated  response  at  the  GaAs  bandgap, 
which  is  independent  of  applied  bias.  The  forward  current  of  units  fabricated 
with  n-type  GaAs  varies  as  exp  (qV/^kT ),  and,  except  at  lower  temperatures 
where  tunneling  becomes  important,  the  values  of  7]  as  a  function  of  the  GaAs 
impurity  doping  concentration  can  be  explained  in  terms  of  the  increase  in  the 
heterojunction  barrier  height  with  voltage,  as  predicted  by  the  graded-gap  model. 
Similarly,  the  voltage  dependence  of  the  reverse  current  is  quantitatively 
explained  by  this  model.  The  results  of  the  capacitance-voltage  measurements 
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are  consistent  with  the  current -voltage  and  photocurrent  measurements 
evaluated  in  terms  of  the  graded-gap  heterojunction  model. 

E.  D.  Hinklcy 
R.  H.  Rcdikcrt 


D.  ION  IMPLANTATION  OF  PHOSPHORUS  IN  GERMANIUM 

A  400  keV  Van  de  Graaf  positive  ion  accelerator  has  been  used  to  implant  phosphorus  ions 
in  germanium.  Experimental  evidence  indicates  that  the  implanted  ions  have  produced  n-type 
regions  in  bulk  p-type  material.  Phosphine  (PH^)  gas  located  at  the  positive  terminal  of  the 
accelerator  passes  through  a  controlled  "leak11  into  a  high-frequency  discharge  tube  where  the 
gas  is  ionized.  The  ionized  components  are  accelerated  along  an  evacuated  column  into  the 
field  of  an  electromagnet  which  analyzes  the  beam  and  deflects  the  phosphorus  ions  onto  the 
germanium  target.  A  polished  [111]  germanium  surface,  inclined  slightly  from  normal  incidence 
in  order  to  reduce  the  probability  of  channeling,  was  exposed  to  the  impinging  phosphorus  beam. 
The  340  keV  phosphorus  beam  was  circular  in  cross-scction  with  a  diameter  of  2  mm.  Implan¬ 
tation  proceeded  for  one  hour  with  an  incident  phosphorus  beam  current  of  0.4  pA,  yielding  a 
1 5 

total  of  9  X  10  implanted  atoms.  Assuming  that  the  phosphorus  atoms  are  distributed  uniformly 

over  a  distance  of  approximately  2  |jl  along  the  direction  of  incidence,  the  concentration  of  piios- 

21  -3 

phorus  atoms  in  the  implanted  region  is  of  the  order  of  1  X  10  cm 

No  observable  macroscopic  surface  damage  was  produced  by  the  implantation,  although 
the  exposed  region  did  become  darker.  Examination  with  an  electron  beam  microprobc  indicated 
that  the  phosphorus  concentration  at  the  sample  surface  was  of  the  order  of  1  X  10^°  cm 
The  implanted  germanium  wafer  was  angle-lapped  and  stained  with  a  CuSO^-HF  solution.  Mi¬ 
croscopic  examination  indicated  that  the  n-type  region  was  approximately  2  p  deep. 

Whereas  the  unimplanted  germanium  showed  p-type  surface  conductivity,  thermal  probe 
measurements  indicated  that  the  implanted  region  was  definitely  n-type.  Post-implantation 
annealing  was  not  required,  although  in  cases  where  it  was  performed  the  n-type  conductivity 
was  enhanced.  A  250 -p  diameter  tin  sphere  was  alloyed  on  the  implanted  side  of  the  germa¬ 
nium  in  the  unexposed  p-type  region  in  order  to  establish  a  "reference1'  p-n  junction,  and  a 
25-(jl  diameter  tin  sphere  was  alloyed  on  the  implanted  n-type  region  to  serve  as  an  ohmic  contact. 
(The  small  size  of  this  latter  sphere  reduced  the  possibility  of  alloying  through  the  narrow  im¬ 
planted  region.)  The  unimplanted  side  of  the  germanium  specimen  was  alloyed  to  an  indium- 
clad  molybdenum  tab  for  ohmic  contact.  Both  p-n  junctions  exhibited  diode  behavior,  but  the 
forward  current  corresponding  to  a  particular  bias  voltage  was  substantially  larger  for  the  ion- 
implanted  junction  than  for  the  "reference"  junction.  This  confirms  that  the  2-mm  diameter 
implanted  region  is  responsible  for  the  observed  germanium  p-n  junction  characteristics. 

E.  D.  Hinklcy 


E.  BULK  GaAs  NEGATIVE  CONDUCTANCE  AMPLIFIERS 

Negative  conductance  microwave  amplifiers  of  high-rcsistivity  n-GaAs  were  first  reported 
by  Thim  ct  al.7  and  have  been  studied  considerably  since. ^  ^  4  Although  it  has  been  assumed 

t  Department  of  Electrical  Engineering  and  Center  for  Materials  Science  and  Engineering,  M.l.T. 
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from  the  outset  that  the  effect  was  due  to  a  bulk  differential  negative  resistance  arising  from 

12  13 

the  transferred-clcctron  mechanism  (which  results  in  Gunn  oscillations  in  lower  resistivity 

material),  this  interpretation  has  never  been  clearly  established.  Computer  calculations  by 

1 1 

McCumber  and  Chynoweth  based  on  their  two-valley  transferred-electron  model  yielded  a 

negative  conductance  at  the  right  frequency,  but  did  not  give  correct  values  of  the  gain  or  oscil- 

14 

lation  threshold.  Similar  inconclusive  results  were  obtained  analytically  using  a  three-slope 
piecewise-lincar  approximation  to  the  drift  velocity-electric  field  relation.  However,  both 
calculations  assumed  a  rather  small  value  for  the  negative  mobility,  whereas  recent  theories  of 
the  electron  transfer^ ^  and  of  the  domain  dynamics  suggest  a  much  larger  negative 

mobility.^ 

By  taking  the  negative  mobility  as  an  adjustable  parameter,  we  have  now  obtained  good 

quantitative  agreement  between  the  measured  small-signal  admittance  parameters  of  GaAs 

negative  conductance  amplifiers  and  those  calculated  from  the  three-slope  piecewise-linear 

2 

model  for  v  vs  E.  The  calculation  yields  a  negative  mobility  of  about  2500  to  3  000  cm  /V  sec, 

which  is  quite  consistent  with  the  transferred-electron  model.  Because  of  the  sensitivity  of  the 

fit  to  the  value  of  the  negative  mobility,  the  above  result  should  be  fairly  accurate  and  perhaps 

constitutes  the  best  present  estimate  of  the  negative  mobility. t 

The  devices  were  prepared  by  alloying  ohmic  tin  contacts  to  opposite  faces  of  chips  of 

4-ohm-cm  n-GaAs,  and  mounting  the  resulting  sample  in  a  small  varactor  diode  package.  Each 

device  was  then  placed  at  the  end  of  a  50-ohm  coaxial  slotted  line,  cooled  to  reduce  its  carrier 

13  -3 

concentration  into  the  proper  range  for  amplification  (about  10  cm  for  a  sample  length  of 

50  p)  and  biased  on  a  pulse  basis  to  avoid  heating.  The  admittance  of  the  samples  was  measured 

with  the  slotted  line  over  a  range  of  frequencies  centered  around  the  frequency  of  maximum 

reflection  gain.  Figure  1-8  shows  the  data  for  one  of  the  devices.  The  sample,  which  was  pre- 

-3  2 

pared  from  n-GaAs  of  4-ohm-cm  room  temperature  resistivity,  was  50  jjl  long  and  1.75  x  10  cm 

13  -3 

in  area.  At  the  operating  temperature  of  — 40°C  the  carrier  concentration  was  1.6  X  10  cm 

2 

and  the  low-field  mobility  was  6000  cm  /V  sec;  the  bias  was  29  V.  The  values  of  the  adjustable 

2 

parameters  used  above  were  \i ^  =  —2500  cm  /V  sec,  -  0,  and  v^  =  0.4  v  .  The  case  capacitance 
of  0.5  pF  (3  00  pBVc  m^)  was  added  to  the  calculated  capacitance  to  give  the  solid  line  shown  in 
the  lower  part  of  the  figure.  Essentially  identical  results  were  obtained  for  the  other  samples 
from  the  same  crystal. 

In  the  upper  part  of  Fig.  1-8,  the  conductance  per  unit  of  sample  area  is  shown  by  the  open 
data  points  for  the  frequency  range  of  2.5  to  5  GHz.  In  the  lower  part,  the  susceptance  is  plotted 
as  an  effective  capacitance  per  unit  of  sample  area  by  dividing  the  measured  susceptance  by  the 
radian  frequency  and  the  sample  area.  Calculated  conductance  and  capacitance  curves  which 
will  be  discussed  later  arc  also  shown.  As  may  be  seen,  there  is  a  broad  region  of  negative 
conductance,  followed  by  a  sharp  rise  to  positive  conductance.  The  capacitance,  on  the  other 
hand,  remains  positive  over  this  frequency  range,  with  a  peak  near  the  zero  crossing  of  the 
conductance. 


tMahrous  and  Robson^  have  calculated  the  impedance  of  GaAs  amplifiers  using  a  three-slope  approximation  to  the  v  vs  E  curve 
of  Butcher  and  Fawcett,  15  but  they  did  not  attempt  a  comparison  with  experiment. 

tGunn  and  Elliott^O  have  recently  inferred  from  very  short  pulse  measurements  on  GaAs  that  the  negative  mobility  is  only  about 
300  cm2/V  sec.  However,  Copeland,  in  a  private  communication,  has  shown  that  their  observations  can  be  explained  with  a 
negative  mobility  of  the  size  found  here  if  the  low-field  hole  generation  rate  which  he  deduced  from  other  experiments^  is  taken 
into  account.  The  extent  to  which  such  hole  generation  would  modify  the  present  results  has  not  yet  been  determined. 


8 


Section  I 


Fig.  1-8.  Comparison  af  measured  small -signal  conductance  and  capacitance  per  unit 
area  (apen  circles)  with  calculations  (solid  lines)  based  on  v  vs  E  curve  af  Fig.  1-9. 


As  a  partial  check  on  the  accuracy  of  the  measurements,  the  data  of  Fig.  1-8  were  used  to 
compute  the  reflection  gain  expected  from  the  device  in  a  50-ohm  reflection  amplifier.  The 
computed  gain  agreed  with  the  measured  values  within  ±0.5dB  over  the  frequency  range  of  2.5 
to  5  GHz. 

We  have  attempted  to  fit  the  data,  using  the  three-slope  piecewise-linear  approximation, 
to  the  v  vs  E  curve  shown  in  Fig.  1-9-  The  values  of  peak  drift  velocity  v  and  peak  electric 

7  P 

field  Ep  were  taken  as  2  x  10  cm/sec  and  4000  V/cm,  respectively,  based  on  threshold  measure¬ 
ments  for  Gunn  effect  oscillators  at  T  =  — 40°C.  The  value  of  p  which  represents  some  average 

2  * 

positive  mobility,  was  thus  taken  as  5000  cm  /V  sec,  slightly  less  than  the  experimental  low- 
2 

field  value  of  6000  cm  /V  sec.  The  mobilities  p^  and  p^  and  the  valley  drift  velocity  v^  were 

considered  as  adjustable  parameters,  although  the  results  are  sensitive  only  to  p?. 

22  Z 
The  spatially  nonuniform  steady- state  solution  may  be  obtained  to  a  good  approximation 

by  neglecting  diffusion  and  assuming  a  virtual-cathode  (E  =  0)  boundary  condition  at  the  negative 

contact.  For  the  AC  solution  we  assume  a  small-signal  perturbation  varying  sinusoidally  in 

time,  assume  that  the  drift  velocity  at  each  point  instantaneously  follows  the  local  electric 
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Fig.  1-9.  Piecewise-linear  approximation  of  drift 
velocity-electric  field  relation  of  n-GaAs  used 
in  admittance  calculation. 


field,!  and  again  neglcet  diffusion.  Under  these  assumptions  and  approximations,  Poisson's 

1 1 

equation  and  the  continuity  equation  combine  to  yield  a  linear  first-order  differential  equation 
for  the  small-signal  eleetrie  field  E^: 


ev(x)  (OE^/Ox)  +  [iu€  +  n(x)  ep(x)]  E^  =  J(cu) 


(1) 


where  v(x)  =  v[E(x)]  is  the  steady-state  electron  velocity,  n(x)  the  steady-state  electron  density, 
|jl  =  dv/dE  the  mobility,  e  the  GaAs  permittivity,  and  J(c o)  the  total  current  density  (including 
displacement  eurrent)  at  frequency  cu.  Solving  Eq.  (1)  with  E^  =  0  at  the  virtual  cathode  at 
x  -  0  and  integrating  the  electric  field  to  obtain  the  small-signal  voltage  V(co),  we  obtain  for  the 
impedanee  of  a  sample  of  length  L 


Z(w) 


V  (cu  ) 

J(w) 


dx^ 


iw[T(X2)-T(x1)] 

e 


ev(xz) 


n(xl )  —  nQ 

n(x2)  -  nQ 


where  nQ  is  the  equilibrium  electron  density  and 


(2) 


t(x)  =  v  *(x)  [x  +  eE/en(x)]  (3) 

is  the  signal  transit  time. 

The  assumption  of  a  pieeewise-linear  v  vs  E  relation  allows  Eq.  (2)  to  be  integrated  easily, 
although  the  resulting  expression  is  rather  involved.!  However,  it  is  a  simple  matter  to  eval¬ 
uate  the  impedance  numerically  for  a  wide  range  of  parameters.  The  calculated  curves  in 
Fig.  1-8  are  for  p^  =  —2500  em  /V  see,  p^  =  0  and  v^  =  0.4  v  .  The  bias  voltage  Vq  and  the 

equilibrium  carrier  concentration  were  set  equal  to  their  experimental  values,  which  were 

13  -3 

Vq  29  V  and  nQ  =  1.6  X  10  em  .  Good  agreement  is  found  for  the  conductance  over  the  fre¬ 
quency  range  2.5  to  4.2  GHz.  Above  4.2  GHz,  the  data  may  be  somewhat  inaceurate  beeause  the 
admittance  measurements  involved  Smith  chart  rotations  of  nearly  one-half  of  a  wavelength 
from  the  reference  short  circuit.  For  the  capacitance,  the  ease  capacity  of  about  0.5pF  (300 pF/ 
em^)  must  be  added  to  the  calculated  curve;  quite  satisfactory  agreement  is  obtained  when  this 
correction  is  made. 


t  The  intervalley  transfer  rate  estimated  in  Ref.  16  is  much  greater  than  our  maximum  frequency  of  5  GHz. 
tlf  fi 3  =0,  the  expression  simplifies  considerably;  the  result  for  this  case  is  given  in  Ref.  20. 
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The  ranges  over  which  the  adjustable  parameters  could  be  varied  with  substantially  the 

2  2 

same  agreement  were  —3000  <  <—2500  cm  /V  sec,  —  250  <  <  250  cm  /V  sec,  and  0.2  <  v  / 

v  <0.5.  Note  that  was  allowed  to  take  on  negative  values  in  order  to  simulate  more  complex 
negative  mobility  characteristics.  As  indicated  by  the  above  ranges,  the  calculated  admittance 
is  sensitive  only  to  the  (initial)  mobility  p^,  which  therefore  permits  a  rather  good  determination 

of  ►‘z- 


The  same  set  of  parameter  values  that  explain  the  admittance  data  also  provide  both  a  rea¬ 
sonable  fit  to  the  measured  current-voltage  characteristics  and  a  prediction  of  the  correct  fre¬ 
quency  and  threshold  bias  for  oscillation. 

Besides  these  results  in  the  3-  to  4-GHz  range,  we  have  some  preliminary  admittance  data 

for  devices  operating  in  the  10-GHz  range  which  can  apparently  be  fit  with  the  same  velocity- 

electric  field  relation.  Also,  this  model  is  able  to  explain  the  very  rapid  variation  of  gain  with 

23 

temperature  reported  by  Thim. 

In  conclusion,  the  pieccwise-linear  model  can  account  for  the  behavior  of  n-GaAs  negative 

conductance  amplifiers  and  in  addition  provides  a  good  estimate  of  the  (initial)  negative  mobility 

for  GaAs.  .  r 

A.  L.  McW  hortcr 

A.  G .  F oyt 


F.  EPITAXIAL  GaAs  GUNN  EFFECT  OSCILLATORS 

Although  Gunn  effect  oscillators  made  from  epitaxial  GaAs  have  been  reported  by  several 
24-27 

workers,  none  of  these  devices  operated  with  power  conversion  efficiencies  as  great  as 

those  obtained  for  devices  made  on  bulk  material.^' ^  We  wish  to  report  a  substantial  improve¬ 
ment  in  the  efficiency  of  such  oscillators  made  with  epitaxial  material.  Efficiencies  as  high  as 
9  percent  have  been  obtained  for  a  device  delivering  a  peak  output  power  of  1.75  watts  at  9.3  GHz. 

The  epitaxial  films  were  deposited  on  GaAs  substrates  using  the  process  described  pre¬ 
viously.  ^  This  process,  first  discussed  by  F.  R.  Knight  et  al.^  and  more  recently  by  other 
32  33 

workers,  '  operates  as  follows:  AsCl^  and  11^  enter  the  hot  zone  of  the  furnace  and  react 

to  form  As  and  HC1.  The  HC1  and  Ga  react  to  form  GaCl,  and  the  As  and  GaCl  are  transported 

down  the  reaction  tube  to  the  substrate.  At  the  (lower)  substrate  temperature,  the  GaCl  reacts 

with  the  As  to  form  GaAs  and  GaCl^.  The  GaAs  deposits  on  the  substrate  as  an  epitaxial  film. 

We  have  had  to  take  two  precautions  in  addition  to  those  already  mentioned  in  the  literature, 

in  order  to  obtain  reproducible  films  suitable  for  Gunn  effect  oscillators.  First,  a  temperature 

gradient  must  be  maintained  in  the  substrate  zone  in  order  to  obtain  a  reproducible  thickness 

of  the  epitaxial  layer. t  Second,  a  prolonged  soak  of  the  substrates  in  hot  potassium  cyanide  is 

necessary  to  prevent  the  occurrence  of  a  p-type  layer  between  the  substrate  and  the  epitaxial 

film.  This  treatment  presumably  removes  copper  from  the  surface  of  the  substrate. 

Some  of  our  best  results  using  this  process  are  shown  in  Table  1-1.  All  these  evaluations 

were  made  on  layers  deposited  on  chromium-doped  semi-insulating  substrates.  As  may  be 

14  15  -3 

seen,  the  best  layers  have  electron  concentrations  in  the  high  10  to  low  10  cm  range, 
both  at  room  temperature  and  at  77 °K.  The  electron  mobilities  for  these  layers  are  typically 


t  For  a  gradient  of  10°C  per  inch,  typical  deposition  rates  are  about  100  p  per  hour  for  a  1 111}  A  face,  and  about  20  p  per  hour  for  a 
{ 100 1  face. 
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TABLE  1-1 

BEST  EPITAXIAL  GROWTH  RESULTS 

n  (300° K) 

n  (77°  K) 

p  (300°  K) 

H  (77°  K) 

Sample 

,  -3, 

(cm  ) 

(cm  3) 

2 

(cm  /V  sec) 

(cm^A /  sec) 

042-2-2 

1 . 2  X  1015 

1 . 2  X  1015 

8,000 

54,000 

043-4 

14 

6. OX  10 

14 

6.  1  X  10 

8,400 

73,000 

043-1-2 

14 

5.  7  X  10 

14 

5.  8  X  10 

8,200 

77,000 

TABLE  1-2 

SUMMARY  OF  OSCILLATOR  RESULTS 

Epitaxial 

Layer 

Thickness 

(h) 

(nl) 

(cm  2) 

Bias 

Vol  tage 

(V) 

Output 

Frequency 

(GHz) 

Peak 

Output 

Power 

(W) 

Effi  ciency 
(percent) 

30 

12 

6X  10 

42 

2.5 

3.2 

7.5 

20 

4X  1012 

11 

5.0 

1.6 

6.5 

10 

2  X  1013 

15 

8.35 

2.5 

9 

10 

2X  1013 

18 

9.3 

1.75 

9 
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2  2 
8000  cm  /V  sec  at  room  temperature,  and  increase  at  77°K  to  values  as  high  as  77,000  cm  /V  sec. 

Wc  have  made  Gunn  effect  oscillators  from  similar  films  deposited  on  heavily  doped  sub¬ 
strates.  The  devices  were  fabricated  by  pla  ing  the  wafer  on  a  tin  clad  molybdenum  tab  and  by 
alloying  the  tin  to  the  GaAs  in  a  reducing  atmosphere.  A  tin  dot  is  then  alloyed  onto  the  epitaxial 
layer  at  temperatures  between  400  and  500°C  in  the  same  reducing  atmosphere.  The  excess 
GaAs  outside  the  area  of  the  tin  dot  is  etched  away,  and  the  resulting  sample  is  mounted  in  a 
standard  varactor  diode  package  for  further  testing.  The  device  is  then  placed  in  a  tunable 
coaxial  cavity  and  biased  on  a  pulse  basis.  The  RF  output  power  is  measured  with  a  thermistor 
power  meter  and  the  output  frequency  monitored  with  a  spectrum  analyzer  through  a  directional 
coupler.  At  each  bias  setting,  the  cavity  was  adjusted  for  peak  output  power.  Some  of  the  best 
results  obtained  to  date  are  summarized  in  Table  1-2. 

The  epitaxial  layer  thickness  for  these  devices  ranged  from  3  0  to  10p  and  the  output  fre¬ 
quencies  ranged  from  2.5  to  9-3  GHz.  Some  of  these  devices  could  be  biased  to  average  electric 
fields  over  10,000  V/cm,  and  gave  up  to  9  percent  efficiency.  We  also  note  that  the  carrier 

concentration  —  sample  length  products  (nl)  for  these  samples  are  considerably  larger  than  the 
12  2 

value  of  about  2X10  /cm  which  has  been  proposed  as  an  optimum  for  oscillator  efficiency. 

In  conclusion  we  feel  that  this  work  has  shown  that  Gunn  oscillators  with  relatively  high 
efficiencies  can  be  fabricated  using  epitaxially  grown  GaAs  and  that  both  the  preparation  of  the 
epitaxial  material  and  the  method  of  contacting,  as  well  as  careful  design  of  the  resonant  cavity, 
are  important  in  the  high  efficiency  which  these  devices  have  shown. 

A.  G.  Foyt 
C.  M.  Wolfe 
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II.  OPTICAL  TECHNIQUES  AND  DEVICES 


A.  OPTIMUM  HETERODYNE  DETECTION  AT  10.6  jum  IN  PHOTOCONDUCT IV E  Ge:Cu 

1-3 

Using  a  Ge:Cu  photoconductor  as  a  heterodyne  detector  at  10.6  pm,  we  have  observed  a 
minimum  detectable  power  which  is  within  a  factor  of  10  of  the  theoretical  quantum  limit,  h^B. 

In  the  experimental  arrangement  used,  the  radiation  from  a  CO^-N^-He  laser,  emitting  approxi¬ 
mately  10  watts  at  10.6pm,  was  incident  on  a  modified  Michelson  interferometer  (see  Fig.  11-1). 
One  mirror  of  the  conventional  interferometer  was  replaced  by  an  off-center  rotating  aluminum 
wheel  with  a  roughened  surface.  The  diffusely  scattered  radiation  from  the  wheel  provided  a 
Doppler  shifted  signal  which  was  recombined  at  the  beam  splitter  with  the  unshifted  (or  local 
oscillator)  radiation  reflected  from  the  mirror  of  the  other  interferometer  leg. 


ROTATING 

WHEEL 


Fig.  II— 1  •  Experimental  arrangement  of  heterodyne  detection  system. 
Electric  field  vector  lies  perpendicular  to  plane  of  paper. 


Figure  11-1  shows  the  experimental  setup.  A  2.54-cm  focal  length  Irtran  II  lens  inserted  in 
the  signal  beam  focused  the  radiation  to  a  single  point  on  the  rim  of  the  rotating  wheel.  The 
function  of  the  lens  was  to  collect  sufficient  scattered  radiation  to  permit  an  incoherent  (non¬ 
heterodyne)  measurement  of  the  scattered  signal  power  at  the  detector  for  calibration  purposes, 

4 

and  to  insure  spatial  coherence  of  the  scattered  radiation  at  the  detector.  Irises  were  used  to 

maintain  the  angular  alignment  of  the  wave  fronts  of  the  two  beams  to  ^2  mrad,  which  is  well 

within  the  required  angular  tolerance  for  optimum  photomixing  (\/a  ~  5  mrad  for  a  detector 
2 

aperture  a  =  2  mm).  A  Perkin-Elmer  wire-grid  polarizer  insured  that  the  recombined  beams 
had  a  common  linear  polarization. 

With  a  bias  voltage  of  13.5  volts  on  the  Ge:Cu  detector,  its  (incoherent)  low-power  respon- 
sivity  was  found  to  be  0.2  amp/ watt  by  calibration  with  a  black-body  source  of  known  temperature. 
The  heterodyne  signal  obtained  from  the  load  resistor  was  fed  into  a  thermocouple-type  rms 
voltmeter  through  a  controlled-bandwidth  low-noise  preamplifier.  Alternately,  the  load  re¬ 
sistor  output  was  fed  simultaneously  to  an  oscilloscope  and  to  a  spectrum  analyzer. 

The  results  of  a  typical  measurement  are  shown  in  Fig.  11-2.  The  solid  line  is  the  observed 
signal-to-noise  power  ratio  (S/N)powpr  of  the  heterodyne  signal  as  a  function  of  the  signal  beam 
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P#  (waits) 


Fig.  11-2.  Theoretical  and  experimental 
values  af  signal-ta-naise  power  ratio 
as  a  function  af  signal  power. 


Fig.  11-3.  (a)  Multiple-sweep  display  af 

heterodyne  signal.  Loss  af  definition  af 
waveform  in  third  cycle  reflects  finite 
bandwidth  af  heterodyne  signal,  (b)  Sin¬ 
gle  sweep  af  heterodyne  signal  shown  in 
(a),  but  with  a  longer  time  scale.  Mod¬ 
ulation  of  signal  envelope  arises  fram 
random  nature  af  scattering  surface. 


H  h 


10  psec 


(b) 


MW 


H  h 


50  psec 
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radiation  power  (Pg)  for  a  typical  run.  Only  noise  arising  from  the  presence  of  the  L.  O.  beam 

(which  was  the  dominant  contribution  to  the  noise)  is  considered.  Various  values  of  P  were 

s 

obtained  by  inserting  calibrated  CaF^  attenuators  in  the  signal  beam,  while  the  unattenuated 
power  was  measured  by  chopping  the  signal  beam  in  the  absence  of  the  L.  O.  With  a  heterodyne 
signal  centered  at  about  70  kHz,  and  an  amplifier  bandwidth  of  270  kHz,  the  experimentally  ob- 


(defined  as  that  signal  beam  power  for  which  the  hetero- 


served  minimum  detectable  power  Pc 

/  S  -14 

dyne  S/N  is  unity)  is  seen  to  be  3.5  x  10  watts.  This  corresponds,  in  a  1  Hz  bandwidth,  to 

-19 

a  minimum  detectable  power  of  1.3  X  10  watt.  Measurements  were  made  with  an  L.  O.  power 

1  3 

of  1.5  mv\  A  plot  of  the  theoretical  expression  *  (S/N)  =  77P  /ZhvAf.  which  includes  the 

^  '  'power  '  s/  j.  £  * 

effect  of  g-r  (generation-recombination)  noise  from  the  photoconductor  is  also  shown  in 

Fig.  11-2.  With  an  assumed  quantum  efficiency  77  =  l/2,  it  is  seen  to  lie  above  the  experimental 

curve  by  a  factor  of  approximately  2,  but  within  the  limit  of  experimental  accuracy.  Had  noise 

from  sources  other  than  the  L.  O.  been  taken  into  account  in  computing  the  S/N,  the  experimental 

curve  would  be  a  factor  of  3  below  the  theoretical.  Although  the  signal-to-noise  ratio  agrees 

5  6 

closely  with  theory,  the  observed  values  of  signal  and  noise  *  are  higher  than  calculated  on  the 
basis  of  the  measured  responsivity.  An  explanation  for  this  observation  is  being  sought  in  a 
possible  dependence  of  the  photoconductor  gain  on  the  frequency  or  on  the  power  of  the  incident 
radiation. 

Figure  11 -3(a)  shows  a  multiple  sweep  display  of  the  heterodyne  signal  obtained  with  a  signal 

-8 

power  of  i.O  x  10  watt.  The  loss  of  definition  of  the  waveform  in  the  third  cycle  reflects  the 
finite  bandwidth  of  the  heterodyne  signal.  Figure  Il-3(b)  shows  a  single  trace  of  this  signal  for 
a  longer  time  scale.  The  modulation  bandwidth  is  caused  by  statistical  fluctuations  of  the  hetero¬ 
dyne  signal  arising  from  the  moving  diffuse  surface  of  the  wheel.  For  the  experimental  con¬ 
figuration  used,  where  the  roughness  of  the  wheel  (~10fim)  is  comparable  to  the  size  of  the 
focused  spot  on  the  wheel  (~50pm),  the  bandwidth  of  the  noise  modulation,  B,  should  be  approxi¬ 
mately  v/d.  Here,  v  is  the  velocity  at  which  the  illuminated  spot  traverses  the  surface,  and 
d  is  the  diameter  of  the  focused  spot  on  the  wheel.  With  v  -  ru>  and  d  ~  f\/D,  B  is  given  ap¬ 
proximately  by  rwD/FA.  Here  v  is  the  tangential  velocity  of  the  wheel  (157  cm/sec),  u  is  its 
angular  velocity  (10  7rsec  ),  and  r  its  radius  (5.05  cm).  F  represents  the  focal  length  of  the 
lens  (2.54  cm),  while  D  is  the  diameter  of  the  radiation  beam  at  the  output  of  the  laser  (~5  mm). 
Evaluating  this,  we  obtain  a  calculated  value  B  ~  30  kHz  which  is  comparable  to  values  observed 
on  the  spectrum  analyzer.  M  c  Teich 

H. J.  Keyes 
R.  H.  Kingston 


B.  CARBON  DIOXIDE  10.6  MICRON  LASER 

A  sealed  off  Brewster  window  laser  of  a  design  similar  to  the  one  mentioned  in  Solid  State 
Research  Report  (1966:2)  was  operated  to  determine  the  useful  life  of  the  gas  mixture  as  a  func¬ 
tion  of  the  power  output  at  10.6  p.  This  water  jacketed  tube  of  10  mm  bore  and  70  cm  discharge 
length  between  kovar  electrodes,  was  used  in  a  one  meter  half  confocal  cavity  with  output  coupling 
through  multiple  layer  coated  lrtran  11  mirror  which  had  a  22  percent  transmission.  The  tube 
and  1.3  liter  reservoir  were  filled  with  a  mixture  of  2torrs  CO^,  2torrs  N^,  and  7  torrs  He. 
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Fig.  11-4.  Output  power  of  a  sealed-off 
10.6  p  CO2  laser  as  a  function  of  op¬ 
erating  time. 


0  20  40  so  eo  100 

OPERATING  TIME  (hours) 

This  mixture  was  not  accurately  adjusted  nor  were  the  optical  components  chosen  for  maximum 
power  efficiency.  The  maximum  radiation  to  electrical  power  efficiency  measured  previously 
was  7  percent  while  this  tube  was  operated  at  3  percent. 

Figure  II —4  shows  the  output  power  as  a  function  of  operating  time  on  the  laser  tube.  The 
tube  output  power  was  essentially  constant  over  98  hours  which  were  accumulated  over  one  month 
with  most  of  the  data  scatter  shown  due  to  tuning  of  the  laser  mirrors.  This  cavity  configuration 
permitted  tuning  over  several  transverse  modes,  and  the  data  indicated  that  a  reasonable  operat¬ 
ing  life  can  be  expected  from  this  type  of  molecular  laser. 

The  visible  spontaneous  emission  from  the  side  of  the  tube  was  monitored  throughout  the 
operating  time  interval,  and  the  determinations  so  far  indicate  that  CO  angstrom  bands,  N ^ 
second  positive  series  bands,  and  CN  cyanogen  violet  bands  were  present.  The  CO  bands  gen¬ 
erally  were  much  more  intense  during  the  beginning  of  the  test,  but  degraded  to  the  intensity 
level  of  the  and  CN  bands  after  50  hours  of  operation.  Only  a  small  decrease  in  intensity 
was  observed  for  the  and  CN  bands  during  the  operation  of  the  tube. 

Unless  there  is  a  cyclic  process,  the  formation  of  CN  and  CO  will  ultimately  deplete  the 
supply  of  CO^  within  the  tube.  An  understanding  of  the  chemical  processes  taking  place  during 
the  discharge  will  permit  a  realistic  appraisal  of  the  expected  operating  life  of  this  molecular 
laser,  and  may  suggest  possible  means  for  extending  its  life  and  power  efficiency. 

R.  J.  Carbone 


18 


Section  II 


REFERENCES 


1.  A.  E.  Siegman,  S.  E.  Harris,  and  B.  J.  McMurtry,  "Optical  Heterodyning 
and  Optical  Demodulation  at  Microwave  Frequencies,"  Optical  Masers, 
edited  by  J.  Fox  (Polytechnic  Press,  John  Wiley  and  Sons,  Inc.,  New  York, 
1963),  p.  511. 

2.  A.  E.  Siegman,  Appl.  Opt.  5,  1588  (1966);  also  Proc.  IEEE  54,  1350  (1966). 

3.  S.  Jacobs  and  P.  Rabinowitz,  "Optical  Heterodyning  with  a  CW  Gaseous 
Laser,"  Quantum  Electronics  III,  edited  by  P.  Grivet  and  N.  Bloembergen, 
(Columbia  University  Press,  New  York,  1964),  p.  481. 

4.  G.  Gould,  S.  F.  Jacobs,  J.  T.  La  Tourrette,  M.  Newstein,  and  P.  Rabinowitz, 
Appl.  Opt.  3,  648  (1964). 

5.  R.  C.  Jones,  Proc.  IRE  47,  1481  (1959). 

6.  A.  van  der  Ziel,  Fluctuation  Phenomena  in  Semi-Conductors,  (Academic 
Press,  New  York,  1959),  pp.  22,  65. 


19 


III.  MATERIALS  RESEARCH 


A.  INTRODUCTION 

The  crystal  structure  of  Nd^O^e  has  been  determined  by  computer  analysis  of  x-ray  powder 
diffraction  data.  The  atomic  arrangement  is  derived  from  the  K^NiF^  structure,  in  which  the 
rare  earth  atoms  exhibit  an  unusual  9-fold  coordination. 

In  order  to  determine  the  cause  of  the  anomaly  observed  in  the  x-ray  scattering  of  ZnSe, 
neutron  diffraction  measurements  have  been  made  on  this  compound  at  room  temperature  and 
liquid  helium  temperature.  The  results  obtained  so  far  are  inconclusive  because  of  difficulties 
in  evaluating  background  intensity  corrections. 

The  phase  diagram  of  Ag^Te  at  high  pressure  has  been  investigated  by  means  of  electrical 
resistivity  and  x-ray  diffraction  measurements.  Evidence  has  been  obtained  for  three  high- 
pressure  phases,  two  of  which  exist  at  room  temperature. 

Hysteresis  of  the  first-order  BB^  &  B31  magnetic  transition  in  MnAs  has  been  studied  by 
measuring  the  resistivity  and  magnetic  susceptibility  as  functions  of  pressure  and  temperature. 
The  experimental  results,  together  with  earlier  data  for  MnAs^  x^x>  show  that  the  transition 
is  due  to  a  large  exchange  striction  in  the  basal  planes  together  with  a  volume-dependent  Weiss 
molecular  field  and  manganese  atomic  moment. 

The  retrograde  solubility  of  Pb  in  PbS  has  been  investigated  by  means  of  Hall  coefficient 
measurements  on  n-type  samples  prepared  by  annealing  and  quenching  single  crystals  contain¬ 
ing  excess  Pb.  The  carrier  concentration  decreases  with  decreasing  annealing  temperature 
from  2.4  X  1019cm'3  for  912°C  to  1.3  X  1018Cm'3  for  350  °C. 

Measurements  of  the  Hall  coefficient  vs  temperature  for  S-doped  GaSb  have  shown  that 
the  lowest  donor  level  of  S  lies  about  0.075  eV  below  the  lowest  conduction  band  at  atmospheric 
pressure.  The  pressure  dependence  of  resistivity  at  room  temperature  shows  that  this  donor 
level  is  associated  primarily  with  the  conduction  band  minima  at  the  (100)  zone  faces. 

The  magnetoresistance  and  Hall  coefficient  of  high-purity  and  V-doped  Ti^O^  single  crystals 
have  been  measured  at  magnetic  fields  up  to  218  kG.  The  high  magnetoresistance  of  the  high- 
purity  samples  at  4.2°K  is  inconsistent  with  a  conventional  one-band  carrier  transport  mechanism. 
The  p-type  conductivity  and  negative  magnetoresistance  of  V-doped  Ti^O^  suggest  that  part  of 
the  V  is  present  in  the  quadrivalent  state. 

Efficient  pulsed  laser  operation  has  been  obtained  in  Nd-doped  YVO.  single  crystals.  The 

q  +3 

threshold  is  only  about  1  joule.  This  low  threshold  is  due  to  the  reduced  Stark  splitting  of  Nd 
ions  in  YVO^  and  to  substantial  transfer  of  absorbed  energy  from  the  lattice  to  the  Nd  ions. 

A  method  employing  PbTe  and  SnTe  as  standards  has  been  developed  for  the  determination 
of  Pb  and  Sn  in  Pb^  ^Sn^Te  alloys  by  electron  microprobe  analysis.  The  results  are  in  good 
agreement  with  those  obtained  by  wet  chemical  techniques. 

B.  CRYSTAL  STRUCTURE  OF  Nd^Te  FROM  X-RAY  POWDER  DATA 

Kent  and  Eick1  have  reported  the  preparation  of  a  family  of  rare  earth  compounds  with  the 
general  formula  R^O^Te,  where  R  is  the  rare  earth  element.  In  measurements  on  a  compacted 
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TABLE 

lll-l 

OBSERVED  VS 

CALCULATED  X-RAY  INTENSITIES  FOR  Nd^Te 

(a  =  ■ 

4.023  A,  < 

:=  12.78  A) 

HKL 

^obs 

^cal  c 

HKL 

^obs 

^calc 

002 

101 

1.0 

0.9 

1.1 

0.9 

316,  307 

323;  2,  0,  10 

17.0 

17.0 

103 

100.0 

100.0 

0,  0,  12;  228 

1.9 

1.0 

110 

39.8 

38.7 

400;  1,  1,  12 

318,  402 

4.7 

4.4 

112 

1.8 

0.9 

309,  411;  2,  1,  11 

3.5 

3.5 

105,  006,  114 

19.4 

19.4 

1,  0,  13;  327 

211 

116,  204 

0.6 

21.7 

0.3 

20.8 

2,  2,  10;  413 

2,  0,  12;  330 

332 

12.2 

12.2 

107,  213 

36.2 

35.4 

0,  0,  14;  415 

215  ,  206  ,  220 

14.7 

15.0 

406,  334 

8.0 

7.4 

109,  301 

6.0 

7.0 

3,  1,  10;  420 

217,  303;  0,  0,  10 
310,  308,  312 

13.8 

14.4 

3,  0,  11;  329 

1,  1,  14;  336 

2,  1,  13;  424 

11.7 

12.5 

305,  226 

S  1 

5.5 

417 

314;  1,  1,  10 

j  .  i 

2,  2,  12;  408 

2.8 

3.2 

219,  321 

1,0,  11 

8.3 

8.0 

1,  0,  12;  426 

2,  0,  14 

8.6 

10.4 
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powder,  they  found  the  density  of  Nd^O^Te  to  be  5.3g/cc.  This  value,  together  with  their  x-ray 

.  o  o 

diffraction  data,  indicated  a  c/a  ratio  of  about  2,  with  a  =  4.023  A  and  c  =  8.519  A.  However, 
these  parameters  lead  to  unreasonable  values  for  the  bond  lengths  in  the  compound,  and  there¬ 
fore  it  seemed  likely  that  the  unit  cell  was  actually  larger. 

In  recent  pycnometric  measurements  made  at  Michigan  State  University,  a  density  of  7.1  g/cc 
was  obtained  for  Nd?0~Te.  The  increased  value  indicated  a  c/a  ratio  of  about  3,  with  a  4.023  A 
and  c  =  12.78  A.  These  new  dimensions,  together  with  the  fact  that  the  revised  Miller  indices 
suggested  a  body-centered  structure,  led  us  to  consider  the  possibility  of  an  atomic  arrangement 
in  Nd^O^Tv;  derived  from  the  K^NiF^  structure,  in  which  the  rare  earth  atoms  exhibit  an  unusual 
9-fold  coordination.  Such  an  arrangement  has  been  confirmed  by  the  results  of  x-ray  diffraction 
measurements  on  a  powdered  sample  of  Nd^O^Te.  The  relative  intensities  were  measured  with 
a  General  Electric  XRD5  diffractometer  using  Ni-filtcred  Cu  radiation.  The  most  probable 
space  group  is  (14/mmm),  and  the  atomic  positions  are; 


+ 

(0 

0 

0, 

1 

2 

1 

2 

Nd 

0 

0 

z, 

0 

0 

z 

in  4e 

Te 

0 

0 

0 

in  2a 

Ox 

0 

1 

2 

1 

4 ' 

1 

2 

0 

1 

4 

in  4c 

The  structure  was  refined  by  using  a  computer  program,  first  employed  to  analyze  data  for 

LaCoO^  (Ref.  2),  which  adjusts  the  structure  parameters  to  obtain  a  minimum  value  for  the  factor 

R,  defined  as  11 1 1  ,  —1  ,  |/z|l  ,  |.  Application  of  this  program  gave  R  =  0.034  and  the 

ods  calc  caic 

following  values  for  the  parameters:  z  =  0.341,  =  0.18,  -  0.54,  and  Bq^  -  0.79,  where 

the  B  parameters  are  atomic  temperature  factors.  The  observed  and  calculated  intensities  arc 
compared  in  Table  111— 1 ,  and  the  structure  obtained  is  shown  in  Fig.  Ill  —  1 . 


[ 3-SS-S493  [ 
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Fig.  Ill-l.  Crystol  structure  of  Nd^O^Te. 
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The  results  of  the  present  investigation  indicate  that  x-ray  powder  data  may  be  used  success¬ 
fully  for  determining  unknown  crystal  structures,  as  well  as  for  the  usual  indexing  and  evaluation 


of  cell  parameters. 


P.  M.  Raccah 
J.  M.  Longo 


C.  NEUTRON  DIFFRACTION  STUDY  OF  ZnSe 

3 

Raccah,  Arnott,  and  Wold  have  recently  obtained  accurate  relative  x-ray  intensity  data  for 
ZnSe  by  means  of  diffractometer  measurements  on  powdered  samples.  They  found  systematic 
discrepancies  between  their  results  and  the  intensities  calculated  from  theoretical  atomic  scatter¬ 
ing  functions.  These  discrepancies  indicate  that  the  electron  distribution  differs  from  that  given 

4 

by  standard  theory.  As  suggested  by  Kaplan  and  Kleiner,  this  difference  might  result  because 
the  usual  scattering  formulas  neglect  overlap  terms  which  could  be  significant  for  a  broad-band 
semiconductor  like  ZnSe,  in  which  the  electrons  are  quite  delocalized.  However,  incorrect 
treatment  of  vibrational  effects  might  also  contribute  to  the  difference  between  theory  and 
experiment. 

In  order  to  investigate  the  role  of  vibrational  effects,  we  have  undertaken  a  neutron  diffrac¬ 
tion  study  of  ZnSe.  In  a  non-magnctic  material  such  as  ZnSe,  the  neutrons  will  be  scattered 
only  by  the  nuclei,  so  that  the  standard  theory  involves  only  three  potentially  unknown  quantities. 
Since  the  standard  theory  uses  the  same  approximations  to  correct  for  vibrational  effects  in 
both  neutron  and  x-ray  diffraction,  success  in  fitting  the  observed  neutron  intensities  by  op¬ 
timizing  these  three  parameters  would  indicate  that  vibrational  effects  are  not  responsible  for 
the  inadequacy  of  the  x-ray  theory. 

The  M.I.T.  nuclear  reactor  has  been  used  to  obtain  several  complete  neutron  diffraction 
patterns  for  ZnSe  at  room  temperature  and  two  complete  patterns  at  liquid  helium  temperature. 
For  both  temperatures,  analysis  of  the  peak  intensities  showed  a  systematic  departure  from 
theory  at  high  scattering  angles.  This  difference  presumably  occurs  because  partial  overlap  of 
the  high-angle  peaks,  even  at  liquid  helium  temperature  where  vibrational  broadening  should  be 
negligible,  results  in  overestimating  the  background  level.  An  attempt  to  determine  the  back¬ 
ground  by  making  a  complete  run  with  the  sample  holder  empty  was  unsuccessful  because  over 
a  third  of  the  background  originated  in  the  sample.  Therefore  the  data  obtained  so  far  are  in¬ 
conclusive.  An  attempt  is  now  being  made  to  develop  a  procedure  for  determining  the  background 
level  by  analyzing  the  peak  shapes.  ^ 

N.  Menyuk 
P.  M.  Raccah 


D.  POLYMORPHISM  IN  Ag2Te  AT  HIGH  PRESSURES  AND  TEMPERATURES 

5  6 

The  polymorphism  of  Ag^Te  at  atmospheric  pressure  is  well  known.  1  With  increasing 
temperature,  stoichiometric  Ag^Tc  is  transformed  successively  from  the  low-temperature 
monoclinic  phase  into  phases  with  face-centered  cubic  and  body-centered  cubic  symmetry.  In 
the  present  investigation,  the  phase  diagram  of  Ag^Te  at  high  pressures  and  temperatures  is 
being  studied  by  means  of  electrical  resistivity  and  x-ray  diffraction  measurements.  Since 
these  properties  are  sensitive  to  deviations  from  stoichiometry,  the  experiments  are  performed 
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on  samples  cut  from  an  n-type  ingot  which  is  nearly  stoichiometric,  since  it  has  a  carrier  con- 

17  -3 

centration  of  only  1.0  X  10  cm 

The  resistivity  at  hydrostatic  pressures  up  to  16.2  kbar  has  been  measured  between  room 
temperature  and  about  360°C  in  a  piston-cylinder  device  which  uses  Dow  Corning  #702  silicone 
oil  as  the  pressure-transmitting  medium.  At  atmospheric  pressure,  in  agreement  with  the  re- 

7 

suits  of  Bottger  and  Meister  for  a  stoichiometric  sample,  there  is  an  abrupt  increase  in  re¬ 
sistivity  at  148  ±  2°C  due  to  the  monoclinic-f. c. c.  transformation.  Application  of  pressure  in¬ 
creases  the  transition  temperature,  as  expected  from  the  changes  in  volume  and  heat  content 
accompanying  the  transformation.  Table  111-2  lists  the  transition  temperatures  observed  in 
repeated  runs  at  5.2,  10.7,  and  l6.2kbar. 


TABLE  III -2 

PRESSURE  DEPENDENCE  OF  TRANSITION  TEMPERATURES  IN  Ag2Te 

Pressure 

Temperoture  of 

Monocl  inic-to-f.  c.  c. 
Transition* 

(°C) 

Temperature  of  Shorp 
Decreose  in  Resistivity* 

(°C) 

1  otm 

148  ±  2 

(Not  obs.  to  525) 

5.2  kbor 

188  ±  2 

239 

10.7  kbor 

205  ±  1 

260 

16.2  kbar 

220  ±  2 

316  ±  4 

*  For  pressures  ot  which  more  thon  one  meosurement  wos  mode,  the  range 
of  observed  temperotures  is  listed. 

At  temperatures  above  the  monoclinic-f. c. c.  transformation,  the  resistivity  at  atmospheric 
pressure  increases  monotonically  up  to  525° C,  the  highest  temperature  investigated.  (Measure 
ments  above  360 °C  were  made  with  the  sample  in  a  quartz  tube  under  argon  gas.)  At  each  of 
the  elevated  pressures,  however,  there  is  a  sharp  decrease  in  resistivity  at  a  higher  tempera¬ 
ture,  as  listed  in  Table  III  -  2 .  This  decrease  probably  results  from  a  transformation  to  a  new 
high-pressure,  high-temperature  phase  not  observed  at  atmospheric  pressure. 

Evidence  for  two  other  high-pressure  phases  was  obtained  from  measurements  at  room 
temperature  under  quasi-hydrostatic  conditions  in  a  tetrahedral-hinge  apparatus.  In  several 
experiments,  the  resistivity  showed  a  sharp  increase  at  21.5  —  22  kbar,  then  a  sharp  decrease 
of  greater  magnitude  at  24.5  —  25  kbar,  after  which  there  was  no  appreciable  change  up  to  about 
40  kbar,  the  highest  pressure  investigated.  These  observations  suggest  the  presence  of  one 
new  phase  between  22  and  25  kbar,  and  a  second  above  25  kbar. 

The  existence  of  these  phases  has  been  confirmed  by  x-ray  studies  with  a  diamond-anvil 

g 

camera  of  the  type  described  by  Piermarini  and  Weir.  At  an  average  anvil  pressure  of  about 
20  kbar,  the  powder  diagrams  contain  lines  due  to  a  new  phase,  in  addition  to  those  of  the  atmos 
pheric  pressure  monoclinic  phase.  At  an  average  anvil  pressure  of  about  30  kbar,  the  lines  of 
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the  monoclinic  phase  are  no  longer  observed,  and  additional  lines  appear.  Since  it  has  not  been 
possible  to  index  the  entire  set  of  lines  obtained  at  this  pressure,  it  seems  likely  that  more  than 
one  phase  is  present. 

The  results  of  x-ray  measurements  made  with  a  diamond-anvil  camera  are  complicated  by 
the  existence  across  the  anvil  faces  of  pressure  gradients  whose  magnitude  varies  with  sample 
thickness.  The  magnitude  of  these  gradients  in  the  present  camera  was  investigated  by  making 
x-ray  measurements  on  finely  powdered  NaCl.  In  these  experiments,  the  sample  pressure  was 

determined  from  the  measured  lattice  constant  a  by  using  the  compressibility  of  NaCl  calculated 
9  ° 

by  Birch.  An  estimate  for  the  pressure  gradient  was  obtained  from  the  uncertainty  in  a  found 
by  least-squares  analysis,  which  reflects  the  broadness  of  the  diffraction  lines.  The  estimated 
gradients  across  the  portion  of  the  sample  exposed  to  the  x-ray  beam  were  3  and  5  kbar  at  average 
anvil  pressures  of  20  and  30  kbar,  respectively.  In  addition,  there  was  considerable  variation 
from  run  to  run  in  the  sample  pressure  for  a  given  average  anvil  pressure.  Thus  the  sample 
pressure  varied  from  20  to  30  kbar  in  four  runs  at  20  kbar  average  anvil  pressure,  and  from 
24  to  44  kbar  in  four  runs  at  25  kbar  average  anvil  pressure.  These  results  indicate  the  necessity 
of  reducing  the  pressure  gradients  in  the  camera  in  order  to  determine  the  structures  of  phases 
with  narrow  high-pressure  stability  ranges.  ^  ^  Banus 

Mary  C.  Finn 


E.  HIGH-PRESSURE  STUDY  OF  MnAs 

The  effect  of  pressure  in  reducing  the  transition  temperature  T  for  the  first-order  B8  —  B31 

1 0  c  1 

transformation  in  MnAs  has  been  reported  previously.  It  was  found  that  above  a  critical  pres¬ 
sure  of  4.5  kbar  the  B31  phase  was  stable  at  all  temperatures.  It  was  also  observed  that  the 
phase  transformation  was  marked  by  an  increasing  degree  of  hysteresis  as  the  temperature  was 
decreased. 

The  pressure-temperature  hysteresis  curve  for  the  reverse  B31  —  B8^  transformation  has 
now  been  determined.  The  results  are  shown  in  Fig.  111-2,  together  with  those  reported  earlier 
for  the  B8^  —  B31  transition.  For  temperatures  down  to  200°K,  the  hysteresis  data  were  ob¬ 
tained  by  observing  the  discontinuity  in  resistivity  due  to  the  reverse  transformation  which  oc¬ 
curred  as  the  pressure  was  reduced  at  constant  temperature.  The  point  at  atmospheric  pressure 
and  138°K  was  obtained  by  magnetic  susceptibility  measurements.  In  this  case,  a  specimen  was 
converted  to  the  B31  phase  by  applying  5  kbar  and  then  cooled  to  77  °K.  The  pressure  was  re¬ 
leased,  the  specimen  was  transferred  without  warming  up  to  the  cold  stage  (at  4.2°K)  of  a 
vibrating-coil  magnetometer,  and  the  magnetization  was  monitored  as  the  sample  was  allowed 
to  warm  up  to  room  temperature.  Initially  the  susceptibility  was  small  and  decreased  with  in¬ 
creasing  temperature,  but  at  138°K  the  magnetization  increased  abruptly  by  a  factor  of  more 

than  50  to  the  value  characteristic  of  the  B8  phase. 

1  11 

The  results  shown  in  Fig.  Ill  - 2,  together  with  those  obtained  in  an  earlier  study  of 
MnASj  x^x>  establish  the  following: 

(1)  There  is  a  dp/dV  >  0  in  the  temperature  interval  Tt  —  AT  <  T  <  Tt, 
where  Tt  =  130°C  is  the  temperature  of  the  second-order  B31  B8j 
transformation,  AT  «  125°C  and  (jl  is  a  manganese  atomic  moment. 

Since  the  thermal  expansion  coefficient  is  ~  2  X  10_4°C'\  this 


26 


Section  Ml 


Fig.  MI-2.  Hysteresis  of  B8^  ^  B3 1 
transition  in  MnAs. 


implies  a  dp/dV  >  0  in  a  critical  molar  volume  range  Vt  -  AV  <  V  <  Vt, 
where  Vt  is  the  molar  volume  at  Tt  and  AV/V  ~  0.025. 

(2)  A  first-order  B8*  ^  B31  transition  at  Tc  occurs  only  if  the  molar 

volume  at  Tc  falls  within  the  critical  range.  Further,  the  fact  that  the 
low-temperature  phase  is  hexagonal,  with  a  discontinuous  expansion  of 
the  basal  planes  on  cooling  through  Tc,  demonstrates  that  there  is  a 
large,  positive  exchange  striction  in  the  basal  planes  if  V  >  Vt  -  AV  at 
Tc.  This  exchange  striction  has  essentially  disappeared  where 
V  <  Vt  -  AV. 

Bean  and  Rodbell12  have  shown  that  a  first-order  transition  can  occur  at  Tc  if 
Tc  =  To[l  +  /3(V  -  Vo)/Vo]  , 

both  the  coefficient  0  and  the  compressibility  are  large,  and  there  is  a  large  AV  at  Tc  due  to 
exchange  striction.  Since  Tc  is  proportional  to  W p  *  ^ ,  where  W  is  the  Weiss  molecular  field 


and  therefore  required  a  large  dW/dV  >  0.  However, 
<  0  and 

>  2 


and  p*  »  4S(S  +  1)  p^,  it  follows  that 


1  dW  2  dp*  > 
P  *  1 W  dV  p*  dV 

Bean  and  Rodbell  assumed  dp*/dV  =  0, 
analysis  of  available  data  gives  dW/dV 

6  <  0  <  22  for  3  > 
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where  and  are  the  atomic  moments  in  the  low-temperature  H8^  and  the  intermediate 
temperature  B31  phases  respectively.  Interpretation  of  the  first-order  phase  change  appears 
to  require^  a  ~  10. 

A  dp/dV  >  0  requires  a  high-spin  *-  low-spin  transition  in  the  critical  molar  volume  interval, 
and  hence 

<d(e  -  €  _)/dV>  AV  >  0.1  eV 
ex  cf 

within  this  interval  AV.  Here  €  and  are  intra-atomic  exchange  and  crystal-field  splittings, 
respectively,  and  0.1  eV  is  taken  as  the  lower  limit  of  the  d-band  width  for  orbitals  of  e^  sym¬ 
metry.  Since 

A(—  e  e)  »  \  f  .  AV/V  »  0.01  eV 
cf  3  cf 

it  follows  that  the  sharp  transition  requires 

Ac  =  <de  /dV>AV>0.1eV 
ex  ex  x 

14 

Stoner  has  pointed  out  that  there  is  a  maximum  bandwidth,  and  hence  a  maximum  overlap 
f 

integral  A^  for  orbitals  on  neighboring  cations,  that  will  support  spontaneous  band  ferro¬ 
magnetism.  Further,  the  bandwidth  for  bonding  orbitals  is  greater  than  that  for  antibonding 
orbitals,  so  that  A^  <  A^,  where  the  subscripts  refer  to  antibonding  and  bonding  orbitals,  re¬ 
spectively.  Thus  the  two  conditions  for  spontaneous  band  ferromagnetism  are 

f  f 

A,  <  A  and  A  .  <  A  <  A,  , 
be  ah  c  b 


where  high-spin  manganese  requires  the  first  and  low-spin  manganese  implies  the  second.  It 
is  concluded  that  the  unusual  occurrence  of  a  dp/ dV  >  0  over  a  small  volume  interval  AV  mani¬ 
fests  the  transition 


A,  <  A  f  ~  A  ,  <  A  f  <  A, 
be  ab  c  b 


f 

and  demonstrates  that  A^  is  quite  sharply  defined. 


J.  B.  Goodenough 
J.  A.  Kafalas 


F.  RETROGRADE  SOLUBILITY  OF  Pb  IN  n-TYPE  PbS 

The  solubility  of  Pb  in  PbS  as  a  function  of  temperature  has  been  investigated  by  means  of 
Hall  coefficient  (Rxr)  measurements  on  samples  annealed  in  argon  at  temperatures  between  350 
and  912°C  and  then  quenched.  The  samples  were  sawed  from  Bridgman-grown  single  crystals, 
supplied  by  Battelle  Memorial  Institute,  which  contain  excess  Pb  incorporated  during  their 
growth  from  undoped,  nominally  stoichiometric  melts.  For  sufficiently  long  annealing  times 
and  sufficiently  rapid  quench  rates,  the  samples  obtained  consist  of  n-type  PbS  saturated  with 
dissolved  Pb  at  the  annealing  temperature,  together  with  a  dispersed  second  phase  of  electrically 
inactive  Pb  microprecipitates.  The  measured  free  carrier  concentration  (n  =  l/eR_-0lH  is 
proportional  to  the  concentration  of  excess  Pb  remaining  in  solution,  since  the  point  defects 
associated  with  the  presence  of  Pb  in  excess  of  the  stoichiometric  composition  are  shallow 
donors  which  are  completely  ionized  at  77  °K. 
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io3/t  (°k) 


Fig.  1 11-3.  Carrier  concentration  (n)  vs  reciprocal  absolute  anneoling 
temperature  far  n-type  PbS. 

The  experimental  results  are  shown  in  Kig.  Ill  —  3,  where  the  logarithm  of  earrier  eoncentra- 
tion  is  plotted  against  reciprocal  absolute  annealing  temperature.  The  solubility  of  Pb  is  retro¬ 
grade.  since  the  carrier  concentration  decreases  with  decreasing  temperature,  from  2.4  x 

AQ  *  7  jO  _  7  7 

10  cm  for  1 185  °K  (10  /T  =  0.844)  to  1.3  x  10  cm  for  623°K  (10  /T  =  1.605).  The  data 
are  well  represented  by  the  straight  line  shown,  which  corresponds  to  the  expression  n  (1.08  x 
1021)  exp  [-(4.34  x  103)/T|  =  (1.08  x  1021)  exp  [-(0.37  eV)/kT[  cm"3. 

The  two  points  in  Fig.  111-3  designated  as  "equilibrated  with  external  phase"  were  obtained 
for  samples  which  were  initially  converted  to  p-type  by  annealing  in  sulfur  vapor  and  then  sat- 
urated  with  Pb  by  equilibrating  at  895  °K  (10 /T  -  1.117)  with  a  liquid- solid  PbQ  ^  mixture. 

The  data  for  these  samples  are  in  good  agreement  with  those  obtained  for  as-grown  samples 
annealed  in  argon.  This  result  shows  that  the  thermodynamic  properties  of  the  internal  micro¬ 
precipitates  in  the  latter  samples  are  essentially  the  same  as  those  of  the  bulk  Pb-rich  external 
phase.  Therefore  it  can  be  concluded  that  the  straight  line  in  Fig.  111-3  gives  the  carrier  con¬ 
centrations  for  PbS  samples  lying  on  the  Pb-saturated  solidus  line  in  equilibrium  with  Pb-rich 
liquid. 

As  shown  in  Fig.  111-3,  the  maximum  carrier  concentration  observed  by  Nensberg  and 

15  19-3 

Petrov  3  in  recent  annealing  experiments  similar  to  the  present  ones  was  only  1.0  x  10  cm 

It  seems  likely  that  their  quenching  rate  was  not  fast  enough  to  prevent  precipitation  in  samples 

annealed  at  the  highest  temperatures  which  they  investigated.  The  same  explanation  probably 

1 6 

accounts  for  the  fact  that  the  maximum  carrier  concentration  observed  by  Bloem  and  Kroger, 

19  -  3 

who  equilibrated  PbS  samples  with  various  partial  pressures  of  sulfur,  was  only  1.6  x  10  cm 

1 5 

At  lower  annealing  temperatures,  Nensberg  and  Petrov  obtained  carrier  concentrations 
systematically  higher  than  the  present  concentrations.  Possibly  they  did  not  anneal  long  enough 
for  their  samples  to  reach  internal  equilibrium,  although  their  annealing  times  would  have  been 
sufficient  to  equilibrate  the  samples  used  in  the  present  experiments. 
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It  is  generally  believed  that  the  predominant  point  defects  in  Pb-rieh  PbS  are  singly  ionized 
sulfur  vacancies.  If  this  is  the  case,  the  deviation  from  stoichiometry  6  in  PbQ 
is  equal  to  (1.31  x  lo”2^)  n,  and  the  straight  line  in  Fig.  111-3  corresponds  to  6  1.41  x  10-2 

exp  [-(0.37  eV)/kT|.  A .  J .  Strauss 

G.  SULFUR  DONOR  LEVEL  IN  GaSb 

Single  crystals  of  GaSb  containing  S  concentrations  of  1017-1018cm  3  have  been  grown 
from  the  melt  by  the  Czochralski  method.  Preparation  of  S-doped  GaSb  has  not  been  reported 
previously.  Figure  111-4  shows  the  temperature  dependence  of  the  Hall  coefficient  at  atmos- 

1 7 

pheric  pressure  for  a  typical  n-type  S-doped  sample,  together  with  data  reported  previously 
for  a  Se-doped  sample  which  has  about  the  same  Hall  coefficient  at  room  temperature.  \\  ith 
decreasing  temperature,  the  Hall  coefficient  of  the  S-doped  sample  increases  rapidly  due  to 
donor  deionization.  This  shows  that  the  S  donors  have  a  finite  ionization  energy.  In  contrast, 
no  deionization  occurs  for  the  Se-doped  sample,  since  the  lowest  Se  donor  level  is  merged  with 
the  conduction  band.  (The  small  decrease  in  Hall  coefficient  with  decreasing  temperature  ob¬ 
served  for  the  Sc-dopcd  sample  is  due  to  the  transfer  of  electrons  into  the  lowest  conduction 
band  minimum  at  k  -  0  from  the  next  higher  (111)  minima,  which  at  atmospheric  pressure  lie 
only  0.08  —  0.09  cV  above  the  k  -  0  minimum.) 

To  determine  the  ionization  energy  of  S  in  GaSb,  computer  calculations  of  Hall  coefficient 
versus  temperature  were  made  for  a  two-conduction-band,  onc-donor-lc vcl  model,  using  known 
parameters  for  the  k  =  0  and  (111)  conduction  bands,  and  optimizing  the  values  of  donor  con¬ 
centration  (Np)  and  ionization  energy  (E^).  The  acceptor  concentration  was  assumed  to  be 


Fig.  111-4.  Hall  coefficient  (R^)  vs  reciprocal 
temperature  for  n-type  GaSb. 
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17  -3 

1X10  em  ,  a  typieal  value  for  GaSb  erystals  grown  from  stoichiomet rie  melts.  The  cal¬ 
culated  curve  best  fitting  the  experimental  results,  which  is  shown  as  a  dashed  line  in  Fig.  111-4, 

17  -3 

was  obtained  for  N^  =  7  x  10  cm  and  =  0.075  eV  relative  to  the  k  =  0  minimum.  In  an 
attempt  to  improve  the  fit  at  low  temperatures,  where  the  experimental  values  lie  below  the 
calculated  eurve,  calculations  were  also  made  for  the  case  of  two  S  levels,  the  second  one  lying 
closer  to  the  k  =  0  band  than  0.075  eV,  but  no  improvement  resulted. 

The  pressure  dependence  of  room  temperature  resistivity  for  the  S-doped  samples  has  been 
measured  by  Kosieki  and  Paul  of  Harvard  University,  who  originally  suggested  the  preparation 
of  S-dopcd  GaSb.  Their  results,  which  will  be  reported  in  detail  elsewhere,  show  that  the  S 

donor  level  with  En  =  0.075  eV  is  primarily  associated  with  the  (100)  conduction  band,  which 

18  19 

lies  at  an  energy  estimated  to  be  about  0.3  —  0.4  eV  above  the  k  =  0  minimum.  *  Evidence  for 
donor  levels  associated  with  higher-lying  conduction  bands  has  been  reported  for  CdTe,  GaAs, 

and  GaAs  P  alloys,  and  evidence  for  a  Se  level  in  GaSb  associated  with  the  (111)  minima  has 

^  ^  20  21 
been  presented  by  Bate  and  by  Kosieki  and  Paul. 

At  temperatures  below  about  100°K,  striking  non-equilibrium  effects  are  observed  in  re¬ 
sistivity  and  Hall  coefficient  measurements  on  S-doped  GaSb.  For  example,  when  a  sample  is 
cooled  rapidly  from  room  temperature  by  immersing  it  in  liquid  nitrogen,  the  carrier  concen¬ 
tration  immediately  after  eooling  is  much  higher  than  the  equilibrium  value  at  77  °K,  frequently 
by  more  than  an  order  of  magnitude.  The  eoncentration  then  decreases  to  the  equilibrium  value 
with  a  time  constant  whieh  depends  strongly  on  temperature,  deereasing  from  about  25  hours 
at  77 °K  to  about  5  minutes  at  90°K.  Although  no  detailed  explanation  for  the  non-equilibrium 
effects  has  been  formulated,  it  seems  probable  that  the  difficulty  in  transferring  electrons  from 
the  conduction  band  to  the  donor  levels  results  from  the  faet  that  the  lowest  conduction  band 
minimum  is  at  k  =  0,  while  the  donor  levels  are  associated  primarily  with  the  (100)  minima. 

A. J.  Strauss 
G.  W.  lseler 

H.  MAGNETORESISTANCE  AND  HALL  COEFFICIENTS  IN  Ti203 

Single  erystal  boules  of  Ti^O^  were  grown  using  speeial  precautions  to  minimize  their  C 
and  N  impurity  eontent.  Commercial  TiO^  and  TiH^  were  mixed  in  appropriate  amounts,  out- 
gassed  at  10  ^torr,  heated  at  1100°C  in  an  atmosphere  of  gettered  argon,  and  pressed  into 
pellets.  The  material  was  melted  in  a  drawn  molybdenum  erueible  under  gettered  argon,  and 
single  erystals  with  an  oxygen/titanium  ratio  x  =  1.512  were  grown  by  the  Czochralski  technique. 
A  V-doped  single  erystal  was  also  grown  from  a  melt  containing  2  atomie-pereent  V^O^.  The 
total  eationie  impurity  eontent  of  the  boules  was  of  the  order  of  180  atomic  ppm  (exclusive  of 
V),  the  major  impurities  being  Si  and  Fe  at  50  and  70  atomie  ppm,  respectively.  The  N  and 
C  concentrations  were  1800  and  240  atomic  ppm,  respectively. 

The  magnetoresistanee  A p/pQ  and  Hall  coefficient  R  of  these  erystals  were  measured  at 
the  National  Magnet  Laboratory  in  magnetic  fields  H  up  to  218  kG.  As  shown  in  Table  II I  -  3 
and  Fig.  111-5,  the  undoped  Ti^O^  samples  exhibited  anisotropy  in  p  and  Ap/p q.  For  current 
at  a  45°  angle  between  the  e  axis  and  basal  plane  of  the  corundum  structure,  magnetoresistanee 
values  exceeding  3.4  were  attained  at  4.2°K  and  II  =  218  kG,  without  any  indication  of  saturation 
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TABLE  111-3 

RESISTIVITY  AND  HALL  COEFFICIENT  OF  Ti^ 


Direction  of 
Current  Flow 

Resistivity 
(ohm -cm) 

Holl  Coefficient 
(cm'Vcoul ) 

4. 2°K 

77°  K 

300°  K 

4. 2°K 

77°K 

300°K 

Undoped  T^O^ 

(somple  1 15) 

[  c  axis 

17,  22 

4. 9  X  10'2 

- 

+1.9  -+3.0 

- 

- 

45°  to  c  axis 

25,  35 

0.13 

0.085 

+4.8 

- 

+0.031 

||  c  ox  i  s 

68 

0. 13 

0.018 

+0.84  - +2.4 

+0.165 

+0.032 

V-doped  T^O^  (somple  110-111) 

1 1  o  oxis 

-4* 

8.3X  10 

8.5  X  10"4 

9.7  X  10"4 

+0.0080 

+0.0040 

+0.0035 

||  c  oxis 

9.5  X  10-4 

9.7  X  10'4 

- 

+0.0045 

+0.0043 

- 

||  b  oxis 

1.2  X  10"3 

1.2  X  10'3 

1.4X  10“3 

+0.0040 

+0.0045 

+0.0039 

*  in  the  I  imit  of  high  magnetic  field.  As  H  O,  R  +0.  028. 
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H  (kG) 


Fig.  111-5.  Magnetaresistance  at  4.2°K  in  undoped  (sample  115),  with  current  flow 

(I)  and  magnetic  field  (H)  along  specified  crystallographic  directions. 

22 

effects.  This  confirms  earlier  results  which  led  to  the  conclusion  that  the  conduction  prop¬ 
erties  of  Ti^O^  at  4.2 °K  cannot  be  interpreted  in  terms  of  a  conventional  one-band  transport 
mechanism.  At  77  °K,  the  magnetoresistance  was  barely  measurable';  this  suggests  that  a  sig¬ 
nificant  change  in  conduction  properties  occurs  between  4.2  and  77 °K.  The  strong  temperature 
dependence  of  resistivity  is  not  ('worthy. 

The  resistivity  and  magnetoresistance  of  V -doped  Ti^O^  are  shown  in  Table  HI- 3  and 
Fig.  111-6.  In  contrast  to  the  undoped  material,  A p/p()  is  negative;  p  remains  reasonably  iso¬ 
tropic,  almost  independent  of  temperature,  and  quite'  small.  With  the  present  sample  the  nega¬ 
tive  magneto  resistance  effects  at  4.2  °K  are  more  pronounced  than  in  the  V -doped  specimens 
23 

studied  earlier,  which  contained  greater  concentrations  of  N  and  C.  The  anisotropy  effects 

in  the  magnetoresistance  are  also  marked.  At  77  °K,  the  negative  magnetoresistance,  though 

small,  is  still  detectable:  A p/p()  —  3.5  X  10  *  at  H  218  kG. 

The  Hall  coefficient  of  the  undoped  material  was  positive  and  small  in  magnitude;  it  could 

be  determined  only  with  great  diffic  ulty  because  the  true  Hall  voltage  amounted  to  less  than 

0.3  percent  of  the  total  transverse  voltage  across  the  probes.  In  the  V -doped  material,  which 

had  essentially  no  spurious  transverse  voltages,  the  Hall  coefficient  was  easily  measurable; 

R  was  again  positive  and  nearly  independent  of  11.  The  values  obtained  correspond  to  a  carrier 

21  3 

concentration  of  about  1.5  x  10  holes/cm  on  a  one-band  model. 

The  origin  of  the  p-type  behavior  of  the  V -doped  samples  has  not  been  firmly  established; 
this  effect  could  be  due  to  deviations  of  sample  composition  from  stoichiometry.  In  one  sample 
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Fig.  MI-6.  Negative  magnetaresistance  at  4.2°K  in  V -doped  U2O3  (sample  110-111), 
with  current  flaw  (I)  and  magnetic  field  (H)  along  specified  crystallographic  directions. 

(not  used  in  the*  investigations  reported  here)  the  stoichiometry  ratio  x  in  Ti  j „y VyOx  was  found 

to  lie  in  the  range  1.517  to  1.527.  This  lends  credence  to  a  recent  suggestion*^  that  part  of  the 

V  in  the  crystal  may  be  present  in  the  quadrivalent  state,  which  would  require  the  ineorporation 

of  excess  oxygon  into  the  lattice.  The  concentration  of  acceptors  due  to  the  presence  of  excess 

oxygen  in  the  above-mentioned  sample  would  exceed  the  concentration  of  electrons  introduced 

by  the  substitution  of  V  for  Ti,  thereby  rendering  the  material  p-type.  Furthermore,  the 

23 

quadrivalent  V  is  also  associated  with  a  magnetic  moment,  as  outlined  previously  this  is 
required  to  explain  the  negative*  magnetoresistance  effects. 

J.  M.  Honig  R.E.  Fahey 

T.B.  Reed  L.  L.  VanZandt 

I.  UNUSUAL  CRYSTAL- FIELD  ENERGY  LEVELS  AND  EFFICIENT 
LASER  PROPERTIES  OF  YV04:Nd 

Optical  absorption  and  laser  studies  have  been  made  on  single  crystals  of  YVO4  grown 
from  iridium  crucibles  in  an  oxyhydrogen  gas-fired  furnace  by  a  modified  Czochralski  technique. 
This  compound  crystallizes  in  a  tetragonal  space  group  of  the  zircon  (ZrSiO^)  type.  All 

rare  earth  vanadates  have  the  same  structure  and  form  solid  solutions  with  YVO^.  Therefore, 
it  will  be  possible  to  investigate  a  variety  of  cross-pumped  laser  systems,  as  in  the  case  of 
yttrium  aluminum  garnet  (YAG). 


*  The.se  cryslals  wore  supplied  b>  II.  M.  Dess  and  S.  II 


Bolin  of  the  Linde  Division,  Union  Carbide  Corporation. 
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Fig.  111-7.  Absorption  levels  of  Nd  -doped 
YAG  ond  YVO4  between  0.6  and  1.0  p 
measured  at  300°K  for  samples  3  mm  thick. 
Fine  structure  represents  transitions  from  ex¬ 
cited  levels  of  ground  stote.  Numbers  be¬ 
low  obsorptians  represent  the  Stark  split 
levels.  For  YAG:Nd  and  YVO^.Nd,  split¬ 
ting  corresponds  to  tetrogonal  ond  cubic  sym¬ 
metry  respectively. 


At  present,  Nd+^-doped  YAG  is  one  of  the  most  efficient  solid  state  lasers.  Figure  111-7 

+  3  +3 

compares  the  absorption  spectrum  of  Nd  -doped  YAG  and  YVO..  The  Nd  absorptions  are 

4  2  ^ 

labeled  according  to  free  ion,  R-S  coupling.  Koningstein  and  Geusic  have  shown  these  term 

*  +3 

designations  are  appropriate  for  YAG:Nd.  For  this  case  the  Nd  ion  is  situated  at  a  tetragonal 

site.  The  tetragonal  distortion,  although  weak,  is  sufficient  to  completely  remove  the  (J  +  1  /2 ) - 

4  4  2 

fold  degeneracy  and  produces  a  pronounced  splitting  of  the  ^5/2'  anc*  **9/2  ^eve^s 

shown  in  Fig.  111-7. 

In  contrast,  very  little  Stark  splitting  is  observed  for  Nd+^  in  YVO^.  This  reduction  in 
splitting  is  apparent  in  the  absorption  spectrum  of  Fig.III-7  (where  the  resolution  is  greater 
than  10  cm"  ).  It  is  present  throughout  the  visible  spectrum  and  persists  at  low  temperatures. 
The  reduction  is  also  shown  by  the  fluorescence  spectrum  for  YVO^:Nd,  which  contains  only 
half  the  lines  present  in  the  YAG:Nd  spectrum.  Most  of  the  emission  in  YV04:Nd  occurs  at 

o 

1.069  p  in  a  line  whose  halfwidth  is  ~5A  at  77  °K.  Consistent  with  the  reduced  splitting,  the 

4 

F  metastable  level  of  YVO^iNd  has  an  oscillator  strength  larger  and  a  radiative  lifetime 
(33  fisec)  smaller,  by  a  factor  of  about  7,  than  for  YAG:Nd. 


*  As  a  result  of  detailed  crystal-field  calculations  these  authors  have  shown  that  Nd+^  in  Y  AG  is  at  a  site  of  tetragonal  symmetry. 
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Fio.  111-8.  Relative  transmission  of  YAG  and  WO4  from  2000 
ta  5500  A.  Spectra  were  measured  at  300°K  for  samples  3  mm 
thick.  Absorption  levels  af  Nd+3  and  Eu+3  are  a|sa  shown. 


Figure  UI-8  shows  the  relative  transmission  of  YAG  and  two  types  of  YVO^.  "Pure"  YVO^ 

has  a  normal  absorption  edge  and  is  colorless.  Photons  of  wavelength  below  4500  A  cause  it  to 

fluoresce  in  several  lines  between  4800  and  5750 A,  as  previously  reported  by  Brixner  and 
27 

Abramson.  This  emission  is  completely  quenched  in  YVO.rNd  at  300 °K  due  to  energy  transfer 

+  3  4 

from  the  lattice  to  the  Nd  ions.  Only  partial  quenching  occurs  at  77° K.  This  implies  that  the 
energy  transfer  processes  arc'  phonon  assisted.  No  evidence  for  energy  transfer  has  been  re¬ 
ported  for  YAG:Nd. 

o 

A  second  type  of  YVO^  has  a  broad  absorption  peaking  near  4200 A  and  is  yellow.  Samples 

of  each  type  are  being  investigated  by  electron  spin  resonance  in  order  to  detect  the  possible 

+  4  1 

presence  and  site  symmetry  of  V  (3d  )  ions.  The  pulsed  laser  threshold  is  lower  in  the  yellow 

O 

material  than  in  the  colorless.  Therefore  the  additional  absorption  at  3500-5000  A  in  the  yellow 
type  also  transfers  energy  to  the  Nd+^ions. 

Low  pulsed  laser  thresholds  ( —  1  .J )  have  been  obtained  in  YVO  :Nd+\  The  low  threshold  is 

due  to:  (1)  the  reduced  Stark  splitting,  and,  (2)  the  substantial  transfer  of  lattice  energy  to  the 

Nd  ions.  T  „ 

J .  R.  O'  Connor 


J.  ELECTRON  MICROPROBE  ANALYSIS 

A  method  has  been  developed  for  the  determination  of  Pb  and  Sn  in  Pb^_xSn  Te  alloys  by 
electron  microprobe  analysis.  Simultaneous  x-ray  intensity  measurements  of  the  fifth  order 

o  o 

Lev  1  line'  of  Pb(A  1.175  A)  and  the  third  order  Lcv^  line  of  Sn(A  =  3.600  A)  are  made  on  separate 
spectrometers.  Several  areas  of  each  alloy  sample  are  examined.  Intensity  measurements  of 
the  Pb  and  Sn  lines  are  also  made  on  pieces  of  PbTe  and  SnTe,  respectively,  which  are  used 
as  standards. 
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TABLE 

ANALYSIS  OF  Pb 

III -4 

.  Sn  Te  ALLOYS 

1  -x  X 

Microprobe  Anolysis 

Wet  Chemical  Anolysis 

(mole  fraction) 

(mole  froction) 

PbTe 

SnTe 

Total 

PbTe 

SnTe 

Totol 

0.968 

0.031 

0.999 

0.986 

0.017 

1.003 

0.943 

0.062 

1.005 

0.940 

0.057 

0.997 

0.799 

0.207 

1.006 

0.797 

0.206 

1.003 

0.620 

0.378 

0.998 

0.608 

0.385 

1.003 

Calibration  curves  for  Pb  and  Sn  in  the  Pb,  Sn  Te  system  were  calculated  by  using  the 

1-x  x  28  29 

backscattered  electron  correction  of  Duncumb  and  Shields,  the  ionization  correction  of  Nelms, 

and  the  absorption  correction  of  Philibert.'*0  The  mole  fraction  of  PbTe  present  in  each  sample 

is  obtained  from  the  calculated  curve  using  the  ratio  of  the  Pb  line  in  the  sample  to  that  in  PbTe, 

Similarly,  the  mole  fraction  of  SnTe  is  obtained  from  the  ratio  of  the  intensity  of  the  Sn  in  the 

sample  to  that  in  SnTe. 

Four  samples  analyzed  with  the  microprobe  were  also  analyzed  for  Pb  and  Sn  by  wet  chem- 
3 1 

ical  techniques.  The  amounts  of  PbTe  and  SnTe  were  calculated  from  the  analyses  assuming 
each  element  to  be  present  entirely  as  the  telluride.  The  results  obtained  by  both  methods  are 
given  in  Table  III -4.  They  are  in  good  agreement  except  for  the  low  Sn  sample  where  the  error 
in  both  types  of  analysis  is  large.  Mary  c  p,jnn 
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IV.  PHYSICS  OF  SOLIDS 


A.  ELECTRONIC  BAND  STRUCTURE 

1.  Current  Pulse  Modulated  Reflectivity  of  Metals 

The  current  pulse  modulated  reflectivity  technique  reported  in  the  last  Solid  State  Research 
Report  (1966:3,  p.  26),  has  been  applied  to  the  study  of  gold  and  niekel.  The  d-band  metals  are 
of  particular  interest  since,  even  though  there  exist  in  the  literature1  considerable  optical  as 
well  as  photoemission  results,  the  interpretations  of  these  measurements  in  terms  of  band  struc¬ 
ture  is  not  yet  conclusive. 

o 

Films,  of  the  order  of  2000  A  thick,  of  gold  and  nickel  were  evaporated  onto  0.150  mm  thick 
glass  substrates;  the  substrates  were  cemented  to  a  liquid  nitrogen  cooled  heat  sink  with  silicone 
grease  which,  for  improvement  of  its  thermal  conductivity,  was  mixed  with  silver  powder. 
Electrical  leads  were  silver-cemented  to  the  samples.  The  eurrent  pulse  rate  was  about  15  Hz, 
with  a  unit  duty  cycle,  and  the  peak  power  dissipated  in  the  samples  was  about  12  watts/cm  ; 
this  power  dissipation  raised  the  temperature  of  the  film  an  average  of  about  40° K  above  the 
base  temperature  which  was  near  77°K.  The  modulation  of  the  temperature  at  15  Hz  was  of  the 
order  of  a  few  degrees.  The  rest  of  the  system  was  conventional. 

The  results  obtained  for  AR/H  on  the  gold  and  nickel  films  are  shown  in  Figs.lV-1  and  1V-2 

along  with  the  reflectivity  R  of  the  same  films.  The  curves,  particularly  the  results  for  gold, 

show  well  defined  structure  out  to  lOeV,  this  being  the  first  time  modulation  measurements 

have  extended  into  the  vacuum  ultraviolet.  For  gold,  the  results  below  5  eV  are  similar  to  those 
2 

of  Garfinkel  etal. ,  who  used  the  piezoreflectivity  method. 

Although  a  detailed  analysis  and  interpretation  of  the  measurements  is  not  possible  at  this 
time,  a  few  points  concerning  the  results  can  be  made.  For  gold,  one  of  the  most  notable  fea¬ 
tures  of  the  data  is  the  sharp  dispersion-like  shape  around  2.5  eV  which  is  only  0.1  eV  wide  from 
the  minimum  in  the  curve  to  the  maximum.  In  this  region  the  behavior  of  the  dielectric  constant 
and  optical  properties  is  the  result  of  both  free  electron  behavior  and  interband  transitions.  The 
onset  of  interband  transitions  is  what  prevents  the  reflectivity  from  dropping  below  33  percent 
in  this  region.  Since  the  free  electron  contribution  to  the  reflectivity  and  dielectric  constant 

ef  1  -  Wp/lw(w  +  i/rc)]  ,  (1) 

where  is  the  plasma  frequency,  depends  on  the  scattering  time  t^,  any  modulation  of  this 
parameter  is  likely  to  show  in  a  AR/R  measurement.  It  seems  possible  then  that  the  initial 
sharp  structure  around  2.5  eV  may  be  due  not  only  to  the  onset  of  L^  L^  transitions  pre¬ 
viously  reported^  but  also  to  the  modulation  of  t(  .  This  sharp  structure  is  absent  in  nickel 
because  the  onset  of  interband  transitions  in  this  metal  occurs  at  a  much  lower  energy  than  in 
gold;  thus  the  free  electron  effects  which  cause  the  sharp  drop  in  the  reflectivity  in  gold  are 
masked  in  nickel. 
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Fig.  IV  — 1 .  Current  modulated  reflectance  AR/R  and  reflectance 
R  far  gold  films  from  2  ta  10  eV. 


Fig.  IV-2.  Current  modulated  reflectance  AR/R  and  reflectance  R  far  nickel 
films  from  2  ta  10  eV. 
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Section  IV 


All  the  previous  papers^  concerned  with  the  optical  properties  and  band  structure  of  d-band 

metals  have  negleeted  any  effects  of  spin-orbit  interaction.  For  the  atomic  state,  the  d  electron 

4 

spin-orbit  energy  is  of  the  order  of  1.8eV  for  gold  and  about  0.3  eV  for  nickel.  Although  it  has 
not  been  calculated,  a  splitting  of  this  order  could  also  be  expected  for  the  solids  and  should  be 
evident  in  the  AR/R  spectrum. 

Finally,  tt  appears  from  photoemission  studies  on  copper,  silver  and  nickel,^  that  direct 

transitions  in  d  band  metals  may  not  be  the  only  source  of  absorption.  Evidence  suggests  that 

for  many  transitions,  only  the  electronic  density  of  states  is  important  and  that  k  conservation 

1 

is  unimportant.  The  relative  importance  of  direct  and  indirect  transitions  is  thus  in  question 
so  it  appears  that  the  interpretation  of  optical  data  is  not  unambiguous  at  this  time.  Nevertheless, 
it  is  hoped  that  information  and  data  obtained  from  new,  more  sensitive,  experimental  methods, 
such  as  the  differential  optical  measurements,  may  help  the  theorist  in  calculating  the  energy 
bands  of  solids.  The  results  presented  here  using  current  modulation  indicate  that  much  more 
structure  is  present  in  the  optical  properties  of  gold  and  nickel  than  one  would  have  believed 
from  the  relatively  uncomplicated  reflectivity  spectra. 

W .  J.  Scouler 
J.  Feinleib 


2.  Current  Pulse  Modulated  Magnetoreflectivity  in  InSb* 

In  addition  to  using  the  electroreflectivity  and  eurrent  pulse  modulated  reflectivity  techniques 

for  studying  metals  and  semiconductors,  we  have  examined  their  use  for  investigating  magnetic 

phenomena  in  semiconductors.  We  have  found  that  the  current  pulse,  or  thermal  modulation, 

appears  well  suited  for  studying  magneto-oscillations  in  the  differential  reflectivity  near  the 

direct  band  edge  in  InSb  cooled  to  liquid  nitrogen  temperatures.  Under  these  conditions,  with 

a  peak  modulation  power  of  0.5  watt  and  a  magnetic  field  of  100  kG,  magneto-oscillations  arc 

observed  across  the  energy  gap  from  near  0.225eV  up  to  about  0.8  cV  with  excellent  signal-to- 

noise  ratio  and  AR/K  ~  10  These  oscillations  correspond  very  closely  with  the  principal 

oscillations  observed  in  magnetoabsorption,  as  well  as  with  theoretical  predictions.  Further 

optimization  of  the  technique  is  needed  before  we  ean  determine  the  relative  merit  of  this  method 

for  investigating  new  materials.  , 

&  J.  Feinleib 

S.  H.  Groves 
C.  Pidgeon^ 


3.  Magnetoreflection  Experiments  in  Arsenic* 

Oscillations  in  the  magnetoreflectivity  of  single  crystal  arsenic  have  been  observed  at  low 
temperatures,  with  the  magnetic  field  H  1 1  S,  where  S  is  the  Poynting  vector.  These  are  the 
first  experimental  observations  in  arsenic  at  infrared  frequencies,  and  provide  information  on 
certain  energy  bands  lying  within  0.5  eV  of  the  Fermi  energy.  With  H||  trigonal  axis,  oscilla¬ 
tions  have  been  observed  only  in  the  energy  range  0.185  <  fico  <  0.375  eV,  and  they  have  been 


*  I'll  is  work  was  carried  out  using  the  high  field  facilities  of  the  National  Magnet  Laboratory,  M.I.  I. 
t  National  Magnet  Laboratory. 
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identified  with  interband  transitions  across  an  energy  gap  of  0.175  eV  and  a  reduced  effective 
cyclotron  mass  of  0.015  mQ.  These  oscillations  are  of  exceptionally  large  amplitude,  even  larger 
than  the  corresponding  effect  in  bismuth.  Nonparabolic  behavior  is  found  for  the  energy  bands 
associated  with  this  interband  transition. 

With  H  binary  axis,  smaller  amplitude  oscillations  are  observed  in  the  energy  range 
0.366  <  ho?  <  0.7  eV.  These  oscillations  are  identified  with  interband  transitions  across  a  band 
gap  of  0.346  eV.  Since  these  oscillations  vanish  for  the  polarization  E||  bisectrix  axis,  the 
energy  band  extrema  are  expected  to  be  located  at  a  high  symmetry  point  in  the  Brillouin  zone. 
The  relation  of  these  results  to  recent  DHVA  measurements,^  and  to  the  Lin-Falicov  band  model' 
is  being  investigated.  M.  S.  Maltz* 

S.  Fischlert 
M.  S.  Dresselhaus 


4.  Optical  Studies  of  Sulfur-Doped  Silicon 

The  following  is  an  abstract  of  a  paper  which  was  presented  at  the  International  Conference 
on  the  Physics  of  Semi-conductors  at  Kyoto,  Japan,  in  September,  1966,  and  which  is  to  be 
published  in  the  Journal  of  the  Physical  Society  of  Japan,  vol.  21  (1966).  This  work  is  a  continua¬ 
tion  of  that  last  reported  on  page  55  of  Solid  State  Research  Report  (1965:2). 

Sulfur,  a  Group  VI  element,  is  a  mono-  or  divalent  donor  in  silicon.  We 
find  that  it  gives  rise  to  at  least  four  different  impurity  centers,  with  ionization 
energies  ranging  from  0.1  to  0.6  eV.  Two  of  these  centers  (at  low  tempera- 

g 

tures)  are  singly  ionized  and  two  are  neutral.  The  spacings  and  the  polar¬ 
ization  dependence  of  the  intensities  arc  wrell  described  in  the  simple  one -valley 
effective  mass  approximation  for  transitions  which  are  hydrogenic  p-like  in 
this  approximation.  The  singly  ionized  centers  have  an  effective  nuclear 
charge  of  +2,  and  therefore  the  spacings  of  the  transitions  to  the  p-like  levels 

9 

are  four  times  as  large  as  for  the  neutral  centers.  There  are,  however, 
observed  features  not  predicted  in  the  simple  effective  mass  approximation, 
and  some  of  these  will  be  discussed.  One  is  a  splitting  of  the  2p()-like  levels 
in  the  ionized  centers.  In  addition  to  the  p-like  transitions,  lines  have  also 
been  observed  which  we  attribute  to  transitions  from  the  ground  state  to  the 
many- valley  orbital  triplet  levels  (s-like  transitions). 

The  study  of  these  centers  under  a  calibrated  uniaxial  stress  gives  much 

information.  The  behavior  of  the  p-like  transitions  for  all  4  centers  gives 

1 0 

directly  the  same  value  of  the  pure  shear  deformation  potential  of  the  con¬ 
duction  band,  E2  =  7.9  eV.  The  behavior  of  the  s-like  transitions  with  stress 
is  characterized  by  deformation  potentials  which  are  much  smaller  than  the 
band  edge  value. 


*  Electrical  Engineering  Department,  M.I.T. 
t  Present  address:  NASA,  Cambridge. 
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The  stress  dependence  helps  eonsiderably  in  the  identification  of  many  of 
the  transitions.  For  three  of  the  eenters,  splittings  are  observed  which  are 
not  due  to  shifts  of  the  conduction  band  minima,  but  whieh  appear  to  be  a  re¬ 
sult  of  the  non-tetrahedral  symmetry  of  the  eenter.  It  is  expeeted  that  the 
uniaxial  stress  measurements  will  enable  the  site  symmetry  of  all  four 
eenters  to  be  identified. 

One  of  the  centers  does  appear  to  have  a  tetrahedral  symmetry.  For  this 
center  we  have  observed  a  pair  of  very  sharp  absorption  lines,  separated  by 
approximately  0.3  meV,  whieh  are  believed  to  represent  transitions  to  the 
orbital  triplet  level  split  by  spin-orbit  coupling.  The  shift  and  splitting  of  the 
transitions  with  stress  and  the  absorption  intensities  with  polarized  light  are 
in  reasonable  agreement  with  such  an  interpretation. 

W.  E.  Krag  H.J.  Zeiger 

W.H.  Kleiner  S.  Fisehler* 

5.  Observation  of  Polaron  Effects  on  the  Intraband  Magnetoabsorption  of  InSb 

In  the  continuing  study  of  polaron  effeets,  it  has  been  found  theoretically  that  the  polaron 

1 1 

cyelotron  energy,  ho;f  ,  is  a  discontinuous  function  of  magnetie  field  near  the  LO  phonon 
energy,  hco^,  at  least  when  a,  the  electron- LO  phonon  coupling  constant,  is  small.  Above 

is  predieted  to  jump  in  energy  by  an  amount  »afia;o  when  compared  to  the  value  obtained 

by  extrapolation  of  the  low-field  plot  of  fio; In  addition,  for  ~  the  oceurrenee  of  an 

eleetron-lattice  resonance  condition  shifts  the  n  -  1  Landau  level  strongly,  as  predieted  from 

1 Z 

theory  and  observed  previously  in  the  interband  magnetoabsorption.  We  have  examined  the 

intraband  magnetoabsorption  of  n-type  InSb  for  photon  energies  from  4  to  26  meV.  Our  results 

1 Z 

clearly  show  the  offset  a  fit ij  as  well  as  the  beginning  of  the  large  effects  reported  earlier  for 
fio;^  «  h (jj o •  We  also  observe  absorption  lines  associated  with  donor  states,  which  show  polaron 
effects.  This  work  constitutes  additional  experimental  confirmation  of  predieted  polaron  effects 
in  InSb  and  provides  an  estimate  of  the  coupling  constant  a. 

D. H.  Dickey 

E. J.  Johnson 
D.  M.  Larsen 

6.  Pressure  Experiments  on  HgSe,  HgTe  and  HgSexTe^-x  Alloys^ 

Experimental  results  on  the  pressure  variation  of  the  thermoeleetrie  power,  a,  for  n-type 
\  1 3 

HgSe,  HgTe  aind  HgSe^Te^  ^  alloys  have  been  reported  which  show  that  with  increasing  pres¬ 
sure  |  cv  |  increases  for  HgSe  and  decreases  for  HgTe.  If  the  samples  studied  were  indeed  ex- 
trinsie  these  results  give  information  about  the  energy  level  ordering  at  k  =  0,  in  particular  that 


*  Present  address:  NASA,  Cambridge. 

tThis  experimental  work  was  performed  at  Harvard  University  in  collaboration  with  Professor  \L  Paul  and  Dr.  W.  M.  DeMeis. 
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E{  r6)  >  K(rg)  for  HgSe  and  E(F^)  >  EfF^)  for  IlgTe.  As  stated  in  the  earlier  report,1  *  this  re¬ 
sult  for  IlgSe  was  unexpected. 

At  this  point  in  this  incomplete  investigation,  two  new  results  of  pressure  experiments  should 
be  reported.  From  the  first  result  it  is  concluded  that  experimental  data  are  needed  to  establish 
whether  or  not  the  samples  previously  used  were  extrinsic.  From  the  second  result,  which  is 
obtained  from  an  optical  pressure  experiment,  it  is  suggested  that  the  ordering  E{Tg)  >  E(  ) 
is  appropriate  for  IlgSe;  this  conclusion  is  opposite'  to  that  obtained  from  thermoelectric  power 
measurements. 

In  the  first  of  these  experiments,  thermoelectric  power  vs  pressure  measurements  were 

17  18  —  3 

made  on  IlgSe  samples  with  electron  concentrations  between  5x10  and  3X10  cm  .  Since 

the  last  report  a  measurement  of  the  Hall  coefficient  vs  pressure  for  a  sample  of  IlgSe  with 
17  3 

5X10  electrons/cm  has  shown  a  change  of  10  percent  in  6  kilobars.  For  this  concentration, 

then,  the  extrinsic  condition,  which  would  give  a  zero  Hall  coefficient  change  with  pressure, 

14 

does  not  hold.  Piotrzkowski  et  al.  claim  that  the  Hall  coetficient  is  independent  of  pressure 

18  3 

for  their  HgTe,  with  1  x  10  electrons/cm  ,  which  shows  a  decrease  of  |  cv  |  with  pressure. 

This  is  about  the  heaviest  doping  of  IlgSe  for  which  a  definite  change  of  <y  with  pressure  can  be 

detected.  Thus,  a  comparison  of  the  Hall  coefficient  vs  pressure  behavior  for  IlgSe  and  HgTe 

18  -  3 

with  1  x  10  cm  carriers  is  at  least  necessary  before  there  can  be  confidence  in  the  simple 
interpretation  of  the  thermo-electric  power  vs  pressure  results. 

The  second  pressure  experiment  has  been  performed  by  Dr.  W.M.  DeMeis  of  Harvard.  He 
has  measured  the  pressure  dependence  of  the  absorption  edge  in  heavily  doped  n-type  IlgSe  and 
HgTe  (the  absorption  edge  corresponds  to  transitions  from  the  top  of  the  valence  band  to  the 
Fermi  level  in  the  conduction  band).  The  results  for  both  materials  can  be  explained  by  assum¬ 
ing  a  pressure  coefficient  of  12  x  10  eV/bar  for  the  F^  -  energy  gap  and  the  energy  ordering 
E(l'g)  >  E(F^).  The  pressure  coefficient  of  this  gap  is  empirically  known  to  fall  in  the  range 
9  to  15  x  10-6  eV/bar  for  the  diamond  and  zinc  blende  semiconductors.  If,  on  the  other  hand, 
ngSe  has  the  ordering  E(  r^)  >  E(Tg)  and  a  valance  band  maximum  not  at  k  0,  then  a  pressure 
coefficient  about  half  of  that  given  above  is  needed  to  obtain  a  zero  or  negative  thermal  energy 
gap;  this  falls  outside  the  usual  range  of  coefficients  for  this  gap.  Although  the  results  of  this 
experiment,  then,  do  not  conclusively  establish  the  ordering  of  levels  at  k  -  0  in  HgSe,  they 
add  to  the  list  of  indirect  evidence  which  favors  the  E(Fg)  >  E(r^)  band  structure. 

S.  H.  Groves 

7.  Analysis  of  Shubnikov-de  Haas  Results  for  HgSe 
1 5 

Whitsett  has  made  a  thorough  experimental  study  of  the  Shubnikov-de  Haas  (SdH)  effect 

in  n-type  HgSe.  An  important  first  result  of  his  work  is  that  the  stoichiometry  of  IlgSe  can  be 

17  18  3 

varied  to  give  1  x  10  to  5x10  homogeneously  distributed  electrons  per  em  .  Thus,  it  is 

probably  the  best  n-type  zinc  blende  material  in  which  to  study  the  SdH  effect. 

Secondly,  by  measuring  the  temperature  dependence  of  the  oscillation  amplitudes,  Whitsett 

found  the  effective  mass  at  the  Fermi  level  for  various  concentrations  and  was  able  to  assign 

values  to  the  two  parameters  which  are  most  important  for  the  conduction  band  eurvature: 

|E(Tg)  —  E(Tg)|  and  the  momentum  matrix  element  P  between  these  states. 


44 


Section  IV 


Thirdly,  Whitsett  showed  that  the  SdH  oscillation  period  was  only  weakly  dependent  upon 
crystal  orientation  with  respect  to  magnetic  field  direction,  as  expected  for  a  conduction  band 
located  at  k  =  0.  However,  a  beat  frequency  was  observed  in  the  higher  concentration  samples, 
and  for  the  magnetic  field  in  the  (110)  plane,  the  beating  was  found  for  all  directions  except  for 
H||  <110)>.  A  nearly  spherical  Fermi  surface  with  bulges  in  the  <lll)>  directions  was  proposed 
by  Whitsett  to  qualitatively  explain  the  beating  behavior.  For  the  magnetic  field  in  the  <11 1> 
or  <100^>  directions,  for  example,  this  Fermi  surface  has  a  maximum  area  off  k^  =  0  and  a 
local  minimum  in  the  area  at  k^  =  0,  where  k.^  is  parallel  to  the  direction  of  the  magnetic  field. 
These  two  external  areas  give  rise  to  two  SdH  frequencies  which  add  together  in  the  conductivity 
to  give  an  average  frequency  term  modulated  by  a  beat  frequency  term  [sin  A  +  sin  R  =  2  sin  \ 

(A  +  B)  cos  \  (A  —  B)].  For  the  magnetic  field  in  the  <110)>  directions,  on  the  other  hand,  the 
cross  sectional  area  has  only  a  maximum  at  k^  =  0  and  a  minimum  of  zero  area  when  the  k^ 
component  of  the  energy  equals  the  Fermi  energy,  thus  giving  rise  to  only  one  SdH  frequency. 

Although  Whitsett' s  model  qualitatively  explains  the  SdH  behavior  of  HgSe,  apparently  no 
attempt  has  been  made  until  now  to  compare  it  with  the  Fermi  surface  constructed  from  the 
theoretical  E  vs  k  expressions  for  small  gap  zinc  blende  materials  such  as  those  derived  by 
Kanc^  £Qr  jnsb.  The  expression  found  by  Kane  for  either  a  or  F^  conduction  band  is  the 
following 

2  2 

E*  =  E'  +  +  a2A'k2  +  b2  [Mk2  +  (L  -  M  -  N)  f.(k)] 

o 

+  C2  [L'k2  -  2(L-  M  -  N)  f  (k)]  ±  nT2  abBf2(k)  (2) 

where 

f  (k)  =  (k^k^  -I-  k^k^  +  k^k^)/k^ 

1  x  y  xz  yz 

and 

f  (k)  =  (k2(k2k2  +  k2k2  +  k\2)  -  9k  2k  2k  2|1 /2/k  . 

2  lxy  xz  yz  x  y  z  ' 

E'  is  the  conduction  band  eigenvalue  which  results  from  the  diagonalization  of  the  Hamiltonian 
between  the  s-like  (T,)  and  p-like  (  ro  +  T-)  basis  which  includes  the  k  •  p  and  the  largest  spin- 
orbit  interaction  terms,  a,  b  and  c  are  the  normalized  coefficients  which  give  the  admixture 
of  s  and  p  state  basis  in  the  eigenvector  corresponding  to  E'  •  A',  L,  M,  N,  L'  and  B  are 
parameters  which  are  determined  by  the  interaction  between  the  higher  conduction  bands  and  the 
s  and  p  bands. 

It  is  convenient  to  separate  Eq.  (2)  into  three  parts:  The  first,  which  includes  Ef,  is  spher¬ 
ically  symmetric.  It  is  the  predominant  part  of  the  energy  and  gives  the  non-parabolic  shape 

of  the  band  but  does  not  contribute  directly  to  the  beating  effects.  The  second  part  of  the  ex- 

2  2 

pression  adds  warping  to  the  energy  surface  and  can  be  written  (b  —  2c  )  (L  —  M  —  N)  f .  (k).  The 

2  2  1 
combination  (b  —  2c  )  (L  —  M  —  N)  contains  opposing  contributions  from  the  k  •  p  matrix  element 

between  the  p-like  components  of  the  conduction  band  and  the  higher  and  conduction 

bands.  The  F^  band  is  thought  to  be  considerably  higher  in  energy  than  the  and  because 


45 


Section  IV 


of  an  energy  difference  denominator  in  (E  —  M  —  N),  the  interaction  with  T gives  the  predominant 
contribution.  The  sign  and  directional  dependence  of  this  term  are  such  that  they  result  in  the 
deviation  from  spherical  symmetry  that  is  required  by  the  Whitsett  model.  The  last  part  of  the 
energy  expression  consists  of  the  ±B  term  above,  and,  in  the  case  of  a  rg  band,  a  term  linear 
in  k  which  is  not  shown  in  Eq.  (2).  For  the  region  of  k-space  of  interest  to  us  this  linear  k 
term  has  the  same  orientation  dependence  as  the  B  term.  These  contributions  to  the  energy 
arise  from  the  lack  of  inversion  symmetry  in  the  zinc  blende  lattice  and,  consequently,  are  zero 
in  the  diamond  lattice  materials.  They  have  been  given  theoretical  attention  but  as  far  as  we 
know  there  have  been  no  experimental  determinations  of  their  size. 

We  have  used  the  following  procedure  to  calculate  cross  sectional  areas  for  a  sample  of 
given  carrier  concentration:  First,  complete  degeneracy  and  spherical  symmetry  are  assumed 
to  find  the  radius  of  the  Fermi  sphere,  k  ,  from  the  carrier  concentration.  Then  the  spherically 

r 

symmetric  part  of  Eq.  (2)  is  used  to  find  the  Fermi  energy  Ep  by  a  numerical  diagonalization 
of  the  interaction  matrix  for  k  =  kr^,  a  spin-orbit  splitting  of  0.5eV,  and  Whitsett' s  two  param- 
eters  |  E(  rg)  —  E(  F^)  |  =  0. 24  eV  and  P  =  7.2  x  10  cV -  cm.  This  process  also  gives  the  co¬ 
efficients  a,  b,  c  at  k  -  kp.  Then,  the  full  expression  in  Eq.  (2)  with  Ep  on  the  left-hand  side 
is  used  to  find  the  magnitude  of  the  k  vector  k^  which  balances  the  equation  for  k  in  a  given 
direction.  Finally,  the  integration  \  k  ^((p)  d  cp  gives  the  cross  sectional  area  for  <p  in  the 

plane  of  interest. 

The  results  of  the  cross  sectional  area  computation  arc  the  following:  the  quantity  (L  —  M  —  N) 

must  be  at  least  three  times  as  big  in  HgSe  as  in  Gc  and  InSb  for  a  rg  conduction  band,  and  six 

times  as  big  for  a  I  ^  conduction  band,  to  explain  the  beating  by  the  second  part  of  Eq.  (2),  or 

the  Whitsett  Fermi  surface.  Only  when  (L  —  M  —  N)  is  made  this  large  does  the  cross  sectional 

area  off  k^  0,  for  k^  aligned  to  one  of  the  <100>  crystal  areas  for  example,  become  greater 

than  the  area  at  -  0.  Our  knowledge  of  the  position  of  the  I  ^  conduction  band  —  the  quantity 

which  largely  determines  the  size  of  (L  —  M  —  N)  —  comes  from  weak  structure  in  the  high  energy 

i  7 

reflectivity,  which  is  tentatively  identified  with  rg  to  F^  optical  transitions.  From  this  data 
the  values  of  3.05,  3.45  and  5.2eV  have  been  assigned  to  the  rg  —  I  ^  energy  separation  in  Gc, 

InSb  and  HgSe,  respectively.  For  L  —  M  —  N  to  be  large  enough  to  explain  the  beating,  the 
Tg  —  T  15  energy  gap  would  have  to  be  ~  1  eV  for  HgSe  for  a  rg  conduction  band  and  even  less  for 
a  conduction  band.  Therefore,  from  what  is  known  about  the  higher  bands,  the  warping  does 
not  seem  sufficient  to  cause  the  beating.  In  addition,  with  the  warping  large  enough  to  produce 
beating,  the  average  or  high  frequency  Sdll  oscillation  shows  a  much  stronger  dependence  on 
magnetic  field  direction  than  was  observed  by  Whitsett. 

With  this  result,  the  third  part  of  Eq.  (2)  must  be  examined  as  a  potential  source  of  the 
beating.  We  find  that  the  difference  in  areas  corresponding  to  the  zince  blende-splitting  of  a 
Tg  conduction  band  can  give  the  observed  frequencies  for  H||  <11 1>  if  B  is  5  to  10  units  of 

in  magnitude.  This  is  not  unreasonable  when  compared  to  theoretical  estimates.  B  has 
to  be  ten  times  as  large  if  the  conduction  band  is  r^.  For  the  magnetic  field  in  a  <11 1>  direc¬ 
tion  the  Fermi  surfaces  corresponding  to  the  split  bands  do  not  touch  at  k^  =  0.  For  H  1 1  <100)> 
there  are  four  <100>  directions  in  the  k^  =  0  plane  where  the  surfaces  touch  and  for  H||  <110)> 
the  k^  =  0  plane  includes  two  <100>  directions  and  four  <11 1>  directions  where  the  surfaces 
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toueh.  However,  the  Whit  sett  data  require  that  beating  be  present  for  Il|f  <100>  but  not  for 

H  <11  0>  .  It  is  not  yet  clear  what  areas  should  be  used  in  these  directions  or  even  if  this 

semielassical  approach  is  valid  for  calculation  of  the  Sdll  frequencies.  Unfortunately,  the 

quantum-mechanical  calculation  of  the  magnetic  energy  levels  ineluding  all  the  terms  in  Eq.  (2) 

1 8 

is  a  difficult  computational  task.  Bell  and  Rodgers  have  recently  performed  sueh  a  computa¬ 
tion  for  one  direction  of  magnetic  field  in  InSb  in  which  the  magnetic  interaction  Hamiltonian 
matrix  of  infinite  order  is  truneated  to  a  250  x  250  matrix.  We  are  currently  looking  for  a 
simpler  treatment  whieh  ean  explain  qualitatively  the  orientational  dependence  of  the  beating. 

It  is  interesting  that  the  beating  effeets  have  not  been  observed  in  gray  tin  samples  with 

1 9 

concentrations  comparable  to  those  of  HgSe.  Gray  tin  has  the  conduction  band  which  en¬ 
hances  the  warping  part  of  Eq.  (2)  but  has  diamond  symmetry  so  that  there  is  no  zero  field  split¬ 
ting  of  the  two-fold  degeneracy.  Most  other  zinc  blende  materials  in  which  the  SdH  effeet  ean 
be  observed  have  rather  than  TH  conduction  bands.  HgTe  is  an  exception  but  the  highly  doped 
material  wc  have  studied  so  far  is  not  sufficiently  homogeneous  to  give  reproducible  beating 
patterns. 

It  is  a  pleasure  to  acknowledge  P.  Trent  for  suggestions  on  several  of  the  computations  used 

in  this  work.  0  TT 

S.  H.  Groves 

P.  W.  Wyatt* 


8.  Conductivity  Calculations  for  Bismuth 

20 

Using  the  density  matrix  formalism,  we  have  calculated  the  zero  field  conductivity  of 

bismuth  assuming  a)  the  holes  may  be  treated  as  free  partieles  with  an  anisotropie  effective 

21  22 

mass,  b)  the  electrons  may  be  treated  by  the  two-band  model  of  Lax,  Cohen  and  Blount,  e) 

23 

the  velocity  matrix  elements  for  the  electrons  are  those  calculated  by  Wolff  based  on  the  two- 
band  model  of  Lax,  Cohen  and  Blount.  The  reflectivity  and  transmission  eurves  calculated  from 
this  conductivity  are  found  to  agree  quite  well  with  the  eurves  for  transmission  displayed  by 
Boyle  and  Brailsford24  and  Boyle  and  Rodgers2^  and  the  absorptance  curve  of  Nanney 2^  This 
calculation  also  confirms  the  idea  that  the  major  contribution  to  the  "core  dielectric  constant" 
eomes  from  the  electrons. 

This  density  matrix  technique  is  now  being  applied  to  obtain  a  detailed  line  shape  of  the 
magnetorcfleetion  in  bismuth,  using  the  magnetic  energy  levels  of  the  two-band  model  and  the 
field  and  frequency  dependent  velocity  matrix  elements  given  by  the  calculation  of  Wolff. 

M.  S.  Maltzt 
M.S.  Dresselhaus 


9.  Fourier  Expansion  for  Energy  Bands  in  Silicon 

Using  the  Fourier  expansion  scheme  reported  earlier  (Solid  State  Research  Report  1966:2, 
p.42),  the  energy  bands  and  frequency  dependent  dielectric  constant  have  been  calculated  for 


*  Present  Address:  Yale  University, 
t  Electrical  Engineering  Department,  M.l.T. 
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Fig.  IV-3.  Energy  bands  of  silicon  calcu¬ 
lated  by  means  of  Fourier  expansion 
scheme;  energy  band  porometers  have  been 
determined  from  experimentol  volues  of 
energy  gaps  ond  effective  mosses. 


Fig.  IV-4.  Comporison  between  colculoted  ond  experimentol  voriotion 
of  the  imoginory  part  of  the  dielectric  constont  in  silicon. 
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silicon.  The  band  parameters  are  determined  exclusively  from  experiment.  The  calculation 

of  the  dielectric  constant  proved  extremely  useful  in  the  identification  of  observed  experimental 

structure  with  characteristic  features  of  the  electronic  structure.  The  silicon  energy  bands  are 

shown  in  Fig.IV-3  and  the  dielectric  constant  calculated  for  these  bands  is  shown  in  Fig.  IV -4. 

27 

For  comparison,  the  experimental  dielectric  constant  of  Philip  and  Ehrcnreich  is  also  shown 

in  Fig.  IV  -4.  _ 

e  G.  Dresselhaus 

M.  S.  Dresselhaus 


B.  MAGNETISM 


1.  High  Temperature  Expansions  for  the  Classical  Heisenberg  Model. 

I.  Spin  Correlation  Function 

The  following  is  an  abstract  of  a  paper  submitted  for  publication  (see  also  Ref.  28). 

The  diagrammatic  representation  of  the  first  eight  coefficients  in  a  high- 
temperature  series  expansion  of  the  zero-field  spin  correlation  function 
(Sf  •  is  presented  for  any  range  of  exchange  interaction  and  for  any 

lattice.  The  present  calculation  exploits  the  order-of-magnitude  simplifications 

which  occur  in  treating  the  quantum-mechanical  spin  operators  in  the  Heisenberg 

1  /2 

model  as  isotropically- interacting  classical  vectors  of  length  [S(S  +  1)]  .  This 

classical  approximation  —  "the  classical  Heisenberg  model"  —  appears  to  be 
excellent  in  the  critical  region,  T  s  T  ,  for  values  of  the  spin  quantum  number 


S  >  1/2.  The  utility  of  the  high -temperature  series  for  is  discussed: 

Besides  providing  an  efficient  method  of  obtaining  the  high-temperature  ex¬ 
pansions  for  all  the  thermodynamic  functions  together,  it  contains  additional 
physical  information.  Such  applications  of  the  high-temperature  series  ex¬ 
pansion  for  will  be  treated  in  subsequent  papers. 

H .  E.  Stanley 


2.  High  Temperature  Expansions  for  the  Classical  Heisenberg  Model. 

II.  Zero- Field  Susceptibility 

Any  extension  of  the  high-temperature  expansion  of  the  zero-field  susceptibility 


X/XCurie  =  1  + 


£  a^J/kT)* 


i=  1 


(3) 


seems  to  be  highly  impractical  because  of  the  enormous  labor  involved.’  Nevertheless,  there 

31 

exists  a  considerable  need  for  more  terms.  Recently  Stanley  and  Kaplan  (SK)  showed  that 
ordcr-of-magnitude  simplifications  occur  when  one  treats  the  non-commuting  quantum-mechanical 


*  “The  labor  of  calculating  just  one  more  term  in  the  series  is  considerably  greater  than  the  labor  of  calculating  all  the  previous 
terms’”  [P.  J.  Wood,  thesis,  King’s  College,  Newcastle-upon-Tyne,  1958  (unpublished)].  The  series  for  general  spin  S  seems 
limited  to  the  six  terms  calculated  by  Rushbrooke  and  Wood. 29  Recently  nine  terms  (ten  for  the  loose-packed  lattices)  have  been 
obtained  for  S  =  1/2  by  methods  practicable  only  for  S  =  1/2,  1  (Ref.  30). 
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spin  operators  oeeurring  in  the  Heisenberg  Hamiltonian  K  =  —  2..J..  S.  •  S.  as  commuting  vectors 
l/2  13  x3  1  J 

of  length  [S(S  +  1)]  /  ,  suggesting  that  many  more  terms  in  the  series  can  be  obtained  "clas¬ 
sically"  than  "quantum-mechanically."  Moreover,  SK  argued  that  useful  results  could  be  ob¬ 
tained  from  this  classical  Heisenberg  model  or  "infinite  spin  approximation,"  since  errors  arc 
small  and  decrease  rapidly  with  S.  For  example,  the  errors  in  the  extrapolated  estimates  of 
Tc/S(S  +  1)  for  cubic  lattices  are  only  ~4  percent  for  S  =  1  and  ~1  percent  for  S  =  5/2.*  Our 
purpose  here  is  to  present  the  results  of  a  calculation  of  the  first  eight  terms  (seven  terms  for 
close-packed  lattices)  in  the  high-temperature  series,  Eq.  (3),  for  a  Heisenberg  Hamiltonian 

with  nearest-neighbor  interactions  only.^  The  a  {£  ^  8)  are  obtained  directly  from  the  first  eight 
f  2 

coefficients  a  in  the  corresponding  high-temperature  series  expansion  for  the  zero-field  static 
spin  correlation  function*^ 


<S, 


v*  -  2 

i=o 


(4) 


by  means  of  the  relation 


a,  -  [XI  JJ*  Nf1 


Z  - 

fg 


fg 


(5) 


fg  fg  — 

where  X  =  S(S  +  1)  and  a  &  =  X-r  cv.s(d).  The  utility  of  the  diagrammatic  representations  is  that 

fa  -  *  d  1  — 

o^&(d)  depends  only  upon  the  topology  of  the  diagram  d,  and  not  upon  the  angles  between  the 
lines.  Thus  for  a  given  order  £  one  need  calculate  only  the  different  a^{r)  where  t  indexes 
the  topological  type  of  diagrams  d.  Although  there  are  many  different  topological  types  of  dia- 

fg 

grams  with  £  straight  lines  and  one  wavy  line  which  might  potentially  contribute  to  a ^  ,  SK 
showed  that  a  /^( r )  is  zero  for  a  large  majority  of  these  types.  Thus  the  summation  cv  /^(d) 

£  t6  fg  fg  1 

in  Eq.(5)  may  be  replaced  by  L  1 g^&(t)  A(t),  and  the  non-zero  o^&(t)  are  given  (for  £  <  8)  in 
Tables  1  through  IV  of  Ref.  28.  Here  t^  denotes  the  number  of  different  topological  types  of 
diagrams  contributing  to  order  i,  and  A(r)  denotes  the  number  of  occurrences  of  the  diagram 
d  of  topological  type  r,  e.g.,  the  number  of  four-sided  polygons  which  can  be  placed  on  a  square 
lattice  is  N;  the  number  of  hexagons,  3N. 

In  this  fashion  we  calculated  the  a.  for  £  -  1,  2,  .  .  .  ,  8  in  terms  of  basic  "lattiee  constants," 

£  34 

which  in  turn  have  been  tabulated  for  the  more  common  two-  and  three-dimensional  lattices. 

The  "general  lattice"  expressions  are  too  lengthy  to  inelude  here,  so  we  instead  present,  in 

Table  1V-1,  the  a^  for  three  three-dimensional  (fee,  bcc,  and  s.c.)  lattices,  and  for  three  two- 

dimensional  (triangular,  square,  and  honeyeomb)  lattices. 


*  The  errors  in  the  extrapolated  estimates  of  TC/S(S  +  l)  for  the  common  two-dimensional  lattices  are  much  larger,  since  here 

Tc(2)/S(S  +  1)  -  2  -  1/S(S  +  1)  (Ref.  32). 


1"  Eight  terms  for  the  fee,  bcc,  and  s.c.  lattices  have  recently  been  published  by  Dood  and  Rushbrooke.33  Our  general-lattice 
results  agree  with  their  calculation  when  we  specialize  to  the  three  cubic  lattices  they  considered 

t  The  diagrammatic  representation  of  the  first  eight  terms  a ^  in  the  high-temperature  series  expansion  of  the  spin  correlation 
function  <Sf  •  Sg>  p  have  been  obtained  recently. 
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Fig.  IV-5.  The  rotios  of  a,/o^o^  ^  of  successive  terms  in  the  susceptibility 
series  for  the  honeycomb,  squore,  and  triangular  plonor  lottices. 
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It  is  perhaps  not  amiss  to  consider  here  the  effect  of  the  additional  terms  upon  two  hotly- 
debated  physical  questions:*  First,  does  a  phase  transition  exist  in  the  two-dimensional  Heisen¬ 
berg  model  with  nearest-neighbor  ferromagnetic  interactions?  The  classical  Heisenberg  model 
is  of  particular  interest  in  connection  with  this  question  because  it  presents  the  same  "dilemma" 

as  that  which  occurs  quantum-mechanically:  Arguments  plausible  at  low  temperatures  lead  to 

37 

a  magnetization  which  is  infinite  in  one  and  in  two  dimensions,  whereas  the  high-temperature 

expansion  coefficients  strongly  suggest  a  phase  transition  at  a  non-zero  temperature  in  two  as 

35 

well  as  in  three  dimensions  (but  not  in  one).  It  was  recently  pointed  out  that  for  some  two- 
dimensional  lattices  the  ratios  a^/a^  ^  of  successive  terms  in  the  high  temperature  series  ex¬ 
pansion  Eq.(3),  when  plotted  against  l/f ,  seem  to  approach  a  straight  line  for  large  l  in  a 

32 

manner  which  is  as  regular  as  for  three-dimensional  cases.  Hence  the  extrapolation  to  £  -  «*> 
and  the  identification  of  the  intercept  with  a  critical  temperature  can  be  made  as  reliably  in 
two  dimensions  as  in  three.  In  Fig.  IV -5  we  plot  the  ratios  a^/a^a^  ^  against  l/f  for  the  classical 
Heisenberg  model.  The  additional  terms  in  these  classical  cases  strengthen  the  conclusion, 
already  drawn  from  the  first  six  terms,  that  there  probably  is  a  phase  transition  for  the  two- 
dimensional  classical  Heisenberg  model. 

The  second  question  is,  does  the  critical  exponent  y  in  the  assumed  form  of  the  divergence 
of  X>  X  ~  (T  —  T  )  ^  as  T  -  Tc+,  have  the  value  4/3  for  the  fee,  bcc,  and  s.c.  lattices?  Previous 
arguments  for  y  -  4/3  have  been  based  upon  the  assumption  that  the  more  classical  the  model, 
the  more  rapidly  the  terms  of  the  high-temperature  series  settle  into  a  smooth  behavior.  This 
assumption,  coupled  with  the  observation  from  the  first  six  terms  that  y(S  =  «)  =  4/3,  has  led 

to  the  proposal  that  y  -  4/ 3  independent  of  spin  and  lattice.  Recently  it  was  suggested  that  this 
3  6 

proposal  is  false,  there  being  instead  a  slow  but  nevertheless  clear  variation  of  y  with  S.  On 
the  basis  of  the  six  terms  available  for  general  S,  no  variation  with  lattice  within  the  class  of 
fee,  bcc,  and  s.c.  was  detected/  With  the  aid  of  the  additional  terms  now  available  for  S  =  °°, 
it  is  possible  to  detect  a  variation  with  lattice  as  well:  for  the  classical  Heisenberg  model, 
y  -  1.33  ±  0.01  for  the  fee  lattice,  with  y  -  1.355  ±  0.01  for  the  bcc  and  s.c.  lattices.  This  new 
result  is  based  upon  the  comparatively  straightforward  and  simple  "slope  method"  for  determin¬ 
ing  the  value  of  y:  If  x  does  diverge  as  T  ^  T  +  with  a  power  law,  then  for  large  £ ,  a^/a^a^  ^  = 
(Tc/Tm)  ^  +  (y  —  l)/i],  where  is  the  ordering  temperature  predicted  by  molecular  field 
theory.  The  slopes  of  the  straight-line  extrapolations  in  Fig.lV-6  for  the  fee,  bcc,  and  s.c. 


lattices  yield  the  values  of  y  given  above. 


H.  K.  Stanley 


3.  Possible  Phase  Transition  for  Two-Dimensional  Heisenberg  Models 

The  following  is  an  abstract  of  a  paper  which  was  presented  at  the  12th  Annual  Conference 
on  Magnetism  and  Magnetic  Materials,  Washington,  D.  C.  and  which  has  been  accepted  for 
publication/^ 


*  These  two  questions  formed  the  basis  of  two  papers^  36  presented  by  the  author  and  T.  A.  Kaplan  at  the  12th  Annual 
Conference  on  Magnetism  and  Magnetic  Materials,  15  —  18  November  1966,  Washington,  D.C.;  see  also  the  following  section, 
B.3,  of  this  report. 

tin  Ref.  36  it  is  argued  that  for  the  spinel  lattice  with  nonmagnetic  A-sites,  the  theoretical  value  of  y  may  differ  from  4/3 
by  as  much  as  50  percent 
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I// 


Fig.  IV-6.  The  ratios  a^/a^a^  ^  af  successive  terms  in  the  susceptibility 
series  far  the  simple,  body-centered,  and  face-centered  cubic  lattices. 


Recently  the  authors  presented  evidence  suggesting  the  presence  of  a  phase 
transition  for  two-dimensional  Heisenberg  models  with  nearest-neighbor  ferro¬ 
magnetic  interactions  and  S  >  Here  we  further  analyze  the  first  six  coefficients 
of  the  high-temperature  series  for  the  zero-field  susceptibility,  and  incorporate 
the  results  of  recent  calculations  beyond  sixth  order  into  the  arguments  for  a 
phase  transition.  It  is  found  that  the  higher  order  coefficients  support  the  pre¬ 
vious  suggestion  concerning  the  existence  of  a  phase  transition  for  S  >  and 
are  consistent  with  a  zero  critical  temperature  for  the  particular  case  S 
Finally,  two  preliminary  calculations  on  the  nature  of  the  low-temperature 
phase  are  described.  The  first  of  these,  due  to  F.  J.  Dyson,  is  an  intuitive 
generalization  of  spin  wave  theory.  The  other  is  an  approximate  calculation 
for  the  rectangular  net  within  the  classical  Heisenberg  model:  Considering 
this  net  to  consist  of  linear  chains,  the  intra- chain  interactions  are  treated 
exactly,  while  the  inter-chain  interactions  are  treated  via  a  molecular  field 
approximation.  Both  of  these  calculations  support  the  original  suggestion  of 
a  phase  transition  for  two-dimensional  Heisenberg  models. 


H.E.  Stanley 
T.  A.  Kaplan 


54 


Section  IV 


4.  Zero  Sound  and  Plasma  Oscillations 

The  properties  of  spin  sound,  the  propagation  of  waves  of  spin  density,  in  unmagnotized  or 
weakly  magnetized  Fermi  gases  have  been  studied  in  the  time  dependent  Hartree- Fock  approxi¬ 
mation,  The  results  have  been  reported  previously  (Solid  State  Research  Report  1965:3,  p.  49) 
and  will  appear  in  Physical  Review.  Two  results  of  particular  interest  were  the  discovery  of 
a  short  wavelength  cutoff  to  the  spectrum  of  the  oscillations  and  a  demonstration  that  portions 
of  the  spectrum  which  were  Landau  damped,  (that  is,  where  single  fermions  were  able  to  absorb 
spin  waves  in  energy  and  momentum  conserving  processes),  nonetheless  should  produce  ex¬ 
perimentally  observable  resonances,  indeed,  of  no  greater  linewidth  than  portions  not  suffering 
from  Landau  damping. 

As  a  sequel  to  this,  the  properties  of  density  oscillations,  ordinary  zero  sound,  have  been 
studied  in  the  same  approximation.  When  the  fermions  interact  through  Coulombic  repulsive 
forces,  as  do  electrons,  we  recover  the  well  known  plasmon  dispersion  equation.  In  addition, 
we  have  found  an  additional  branch  of  the  density  oscillation  spectrum  in  this  approximation,  a 
set  of  roots  to  the  well  known  dispersion  equation  which  has  been  heretofore  neglected.  The 
oscillations  of  this  second  branch  are  Landau  damped,  but  because  they  occur  at  low  frequency 
and  with  a  phase  velocity  very  near  the  Fermi  velocity,  they,  like  the  Landau  damped  spin 
sound  modes,  could  lead  to  experimentally  observable  resonances.  The  short  wavelength 
cutoff  is  also  characteristic  of  density  oscillations.  In  a  dense  Coulomb  gas,  the  cutoff  occurs 
at  a  frequency  above  the  limits  of  definition  of  the  modes  and  is  hence  only  a  mathematical 
quirk.  A  Fermi  gas  with  short  range  interparticle  forces,  however,  offers  some  possibilities 
for  observing  the  cutoff  at  relatively  lower  frequencies  where  the  entire  spectrum  is  well 
defined.  The  spectrum  has  been  studied  using  parameters  of  the  interactions  appropriate  to 

3 

liquid  He  .  The  results  of  this  work  have  been  submitted  for  publication. 

L.  L.  VanZandt 

C.  TRANSPORT  THEORY 

1.  Transport  Theory  of  Interacting  Electrons 
in  the  Random  Phase  Approximation 

In  order  to  study  the  influence  of  the  electron-electron  interaction  on  the  transport  properties 
of  the  conduction  electrons  in  a  metal,  we  have  considered  the  model  of  an  interacting  electron 
gas  driven  by  an  external  electromagnetic  field  and  scattered  by  randomly  distributed  impurities. 
As  a  first  step  for  a  systematic  many-body  theory  of  this  problem,  we  treat  the  electron-electron 
interaction  in  the  random  phase  approximation,  which  is  valid  for  sufficiently  high  electron 
density.  Thus  this  theory  includes  all  the  plasma  (collective)  effects  on  the  electromagnetic 
properties. 

We  have  considered  both  longitudinal  and  transverse  external  electromagnetic  fields.  We 

express  the  linear  response  of  the  system  in  terms  of  the  average  total  electromagnetic  field, 

by  constructing  a  generalized  kinetic  equation.  To  achieve  this  we  have  made  use  of  the  tech- 

38 

nique  we  developed  earlier  for  the  case  of  non-interacting  electrons.  This  kinetic  equation 
for  the  impurity-averaged  density  matrix  is  valid  for  arbitrary  frequencies  of  the  external  field 
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and  for  a  general  onc-clectron  Hamiltonian.  In  addition,  the  effects  of  the  electron-impurity 
interaction  are  described  without  any  approximation.  The  Born  approximation  for  the  scattering 
is  easily  derived. 

The  effects  of  the  electron-electron  interaction  appear  in  various  ways.  First,  the  impurity 
potential  is  screened  by  the  gas  in  thermal  equilibrium,  i.e.,  statically,  and  the  external  field 
is  screened  by  the  dielectric  constant  at  the  driving  frequency.  Additional  terms  describe  the 
modification  of  the  scattering  and  interference  terms  of  the  transport  equation. 

We  have  so  far  examined  the  consequences  of  this  theory  for  the  special  case  of  the  homo¬ 
geneous  electron  gas  in  the  Born  approximation  for  the  scattering  and  for  a  homogeneous  driving 
field.  For  w  =  0  we  have  found  that  the  simple  minded  picture  of  an  independent  electron  (quasi¬ 
particle)  driven  by  the  average  total  field  and  scattered  by  a  statically  screened  impurity  is 
valid.  For  a;  7^  0  this  is  no  longer  true;  effects  of  the  dynamic  screening  of  the  impurities  begin 
to  come  in.  For  cut  »  1,  where  r  =  relaxation  time,  again  the  model  of  an  independent  electron 
is  valid  but  now  with  the  impurity  potential  screened  by  the  geometric  mean  of  the  dielectric 

constants  at  frequencies  zero  and  w.  This  latter  description  is  a  very  simple  expression  of 

39 

the  results  of  various  theories  for  this  particular  case.  We  should  point  out  that  these  the- 
39 

ories  are  not  only  restricted  to  the  case  of  uniform  fields,  but  are  also  unsatisfactory  in  that 
the  identification  of  the  total  field  is  incomplete.  At  least  for  the  longitudinal  case  they  give  the 
correct  result  through  the  cancellation  of  two  errors. 

P.  N .  Argyres 


2.  Landau  Theory  of  a  Fermi  Liquid  in  Random  Impurities 


The  transport  theory  for  the  Landau  quasi-particle  distribution  function  in  the  presence  of 
random  impurities  is  incomplete.  Strictly  speaking  it  is  valid  only  for  (jj  -  0,  where  w  is  the 
frequency  of  an  external  driving  field,  and  approximately  correct  for  cot  «  1,  where  t  is  a 
measure  of  the  relaxation  time.  For  c or  »  1,  a  region  of  experimental  interest,  it  is  incorrect. 
Furthermore,  the  "memory"  effects  of  the  collision  processes  have  been  ignored.  We  have 
begun  to  examine  the  possibility  of  a  new  approach  to  this  problem  that  would  allow  us  to  derive 
a  transport  equation  for  the  quasi-particle  distribution  function  for  arbitrary  a>,  subject  only 
to  the  usual  restriction  of  the  existence  of  the  quasi-particlc.  The  essence  of  the  approach  is 
to  treat  the  quasi-particle  as  a  quantum-mechanical  particle  and  derive  first  its  "Schrodinger 
equation"  and  then  its  transport  equation.  p  re  g 

J.  Sigel*gyre 


D.  SCATTERING  EXPERIMENTS  WITH  LASERS 

1.  Stimulated  Brillouin  and  Raman  Scattering 
in  Quartz  Between  2.1  and  293 °K 

Further  results  on  stimulated  Brillouin  scattering  experiments  in  quartz  at  2.1  °K  have  lent 

40 

additional  support  to  the  contention  that  the  anomalous  Brillouin  shifts  reported  by  Krivokhizha 
et  al.  were  misinterpreted  and  are  actually  due  to  intense  stimulated  Raman  radiation  which  sets 


*  Physics  Department,  M  l  T. 
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in  at  low  temperatures.  We  had  previously  observed  a  weak  additional  line  in  radiation  scattered 

out  of  the  x-direction  which,  because  of  its  relative  displacement  (3  percent)  with  respect  to 

_  \ 

successive  interference  orders  of  the  etalon,  was  inferred  to  be  the  466  cm  Raman  line  of 
quartz  (see  Fig.  1V-7). 

In  scattering  along  the  z-axis,  we  observe  a  shift  (with  or  without  the  regular  Brillouin  shift 

-1 

appearing)  which  is  compatible  with  the  128  cm  Raman  line  of  quartz  (see  Fig.lV-8).  Since 
interferometer  measurements  cannot  uniquely  identify  widely  shifted  lines,  we  simultaneously 
photographed  the  scattered  radiation  on  a  grating  spectrograph.  In  every  case  the  suspected 
Raman  lines  were  observed  and  showed  complete  correlation  with  the  unidentified  Fabry-Perot 
etalon  lines.  For  example,  the  spectrogram  taken  simultaneously  with  the  interferometer  picture 
of  Fig.  IV -8  shows  only  a  strong  line  at  132  cm  and  no  other  lines.  Thus  the  subsidiary  line 
in  Fig.lV-8  is  identified  as  "128  cm  "  stimulated  vibration  of  the  classical  room  temperature 

Raman  spectrum  (duly  temperature  shifted).  Consequently,  at  2.1  °K,  for  z-directed  scattering, 

-1  - 1 

the  Raman  threshold  for  the  128  cm  line  is  lower  than  that  for  the  466  cm  line;  for  x-directed 

- 1  -1 

scattering  both  the  466  cm  and  128  cm  line  often  appear.  In  all  cases  the  stimulated  Raman 

radiation  is  sufficiently  intense  to  pass  through  the  narrow-band  ruby-pass  filter  in  front  of  the 

etalon.  There  seems  to  be  a  tendency  for  stimulated  Brillouin  scattering  to  be  suppressed  at 

the  expense  of  stimulated  Raman  radiation.  Such  interference  between  Brillouin  and  Raman 

41 

processes  has  been  observed  previously. 

In  order  to  verify  that  additional  nonlinear  effects  did  not  set  in  at  higher  powers,  we  raised 
the  power  incident  on  the  sample  in  small  increments.  The  measured  results  continued  to  be  the 
same  up  to  lOOMw,  at  which  point  catastrophic  failure  of  the  helium  dewar  occurred. 

Thus  it  appears  that  the  anomalous  Brillouin  shifts  reported  in  Ref.  40  can  be  attributed  to 
stimulated  Raman  scattering,  which  becomes  quite  strong  at  low  temperatures."  This  phase  of 
the  experiment  has  been  completed,  and  the  results  have  been  submitted  to  Science. 

P.  E.  Tannenwald 
J.  B.  Thaxter 
F.  H.  Perry 


2.  Temperature  Dependence,  Orientation  Correlation, 

and  Molecular  Fields  in  Second  Harmonic  Light  Scattering 
from  Liquids  and  Gases 

The  following  is  an  abstract  of  a  paper  which  is  being  prepared  for  publication. 

The  properties  of  second  harmonic  light  scattering  in  fluids  are  derived 

in  terms  of  the  molecular  polarizabilities,  orientation  correlation  among 

molecules,  and  the  electric  field  F  in  the  unilluminated  liquid.  Measurements 

a 

of  scattering  intensity  in  CCl^  and  water,  from  about  10°C  to  60°C,  are  re¬ 
ported.  The  variations  with  temperature,  and  therefore  apparently  the  effects 
of  changes  of  orientation  correlation,  are  slight.  The  influence  on  the  scattering 


*  The  appearance  of  Raman  lines  on  the  habry-Perot  interferogram  enables  us  to  measure  Raman  shifts  and  their  temperature  de 
pendences  to  an  accuracy  of  0.1  cm'l 
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Fig.  IV -7.  Fabry-Perat  interferometer  spectrum  af  backward  Brillauin  scattering 
along  x-directian  in  quartz  at  2.1°K.  Brillouin  shift  is  0.84  cm“^.  The  weak 
line  is  displaced  successively  by  0.11  cm“^  with  respect  ta  each  interarder 
spacing  and  is  identified  ta  be  the  468  cm"^  Raman  line. 


Fig.  IV-8.  Backward  Brillauin  scattering  along  z -direction  in  quartz  at  2.  1  °K. 
Na  Brillauin  shift.  The  weak  line  indicates  an  apparent  Brillauin  shift  af 
2.30  cm-^,  but  is  actually  the  132  cm'^  Raman  line. 
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of  preferred  orientation  in  CCl^,  of  the  type  deduced  from  X-ray  diffraction, 
is  calculated.  The  third  rank  molecular  polarizability  tensor  consists  of  the 
fourth  rank  polarizability  contracted  with  F  ,  in  addition  to  the  term  pre¬ 
viously  known,  which  is  independent  of  F  .  With  F^  arising  from  molecular 
multipole  moments,  the  new  term  explains  some  experimental  depolarizations 
in  liquids,  and  accounts  for  most  of  the  line,  but  not  the  background,  observed 
in  methane  gas.  In  the  case  of  molecular  dipole  moments,  it  can  cause 

spectral  narrowing.  ...  .  .  . 

1  D.  L.  Weinberg 


3,  The  Raman  Spectrum  of  GaS  and  GaSe 

o 

We  have  observed  the  Raman  spectrum  of  GaS  and  GaSe  using  6328  A  and  1.06  p  lasers.  The 
first  order  spectrum  for  each  compound  consists  of  five  lines,  which  occur  for  GaS  at  69.6, 
183.9,285.7,  289.3  and  357.6cm"1,  and  for  GaSe  at  59.6,  133.8,  209.5,  253.8  and  308.6  cm'1 . 
These  compounds  are  layer  structures  in  which  the  layers  are  relatively  weakly  bound  to  each 
other.  Each  layer  has  point  group  symmetry  and  contributes  four  atoms  per  unit  cell.  The 
stacking  of  layers  in  GaS,  and  one  modification  of  GaSe,  retains  the  symmetry  and  results 
in  eight  atoms  per  unit  cell,  while  the  rhombohedral  modification  of  GaSe  has  twelve  atoms  per 
unit  cell,  with  C3v  point  group  symmetry.  Q  B  Wright 

A.  Mooradian 


4.  The  Raman  Spectrum  of  Trigonal  Selenium 

The  Raman  spectrum  of  trigonal  selenium  has  been  observed  using  a  1.06  micron  YAG:Nj*  + 

laser.  The  point  group  symmetry  predicts  a  singly  degenerate  A^  and  two  doubly  degenerate 

A Raman  active  optical  modes,  as  well  as  a  Raman  inactive  A?  optical  mode.  We  observe 

^  -1 

three  strong  first  order  lines  which  we  identify  as  an  A^  mode  at  142.9  cm  ,  and  an  A^  and  A^ 

mode  at  233.7  and  237.0  cm  \  respectively.  This  ordering  is  very  similar  to  that  calculated 
42 

by  Ilulin  for  tellurium  using  elastic  constants  and  a  simple  force  model.  Second  order  lines 
are  observed  at  457,  438,  273,  206,  183,  and  104  cm  * 

A.  Mooradian 
G.  B.  W right 

5.  Interstellar  OH  Maser 

A  Physical  Review  Letter  entitled  "Maser  Model  for  Interstellar  OH  Microwave  Emission" 
appeared  in  the  August  10th  issue.  In  this  letter,  wc  describe  a  model  of  ultraviolet  pumping, 
microwave  saturation,  and  nonlinear  interaction  of  propagating  modes.  The  nearly- circular 
polarization  observed  was  attributed  to  nonlinear  mode  suppression.  More  detailed  computer 
calculations  by  P.  H.  Trent  are  almost  complete  for  the  self-consistent  populations  in  the  OH 
hype rfine -split  A-doublet  levels  under  somewhat  directional  UV  pumping,  taking  into  account 
selective  absorption  in  the  OH  cloud.  Also,  the  effect  of  microwave  saturation  on  mode  sup¬ 
pression  is  being  further  explored,  both  analytically  and  on  the  TX-2  computer. 

M.  M.  Litvak  M.  L.  Meeks  (Group  31) 

A.  L.  McWhorter  H.J.  Zeiger 
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